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Abstract 



o 
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Qh, We report Suzaku observations of the northern half of the Hydra A cluster out to ~ 1 .4 Mpc, reaching the virial 

radius. This is the first Suzaku observations of a medium-size {kT ~3 keV) cluster out to the virial radius. Two 
^ ' observations were conducted, north-west and north-east offsets, which continue in a filament direction and a void 

jXi ' direction of the large-scale structure of the Universe, respectively. The X-ray emission and distribution of galaxies 

C^ \ elongate in the filament direction. The temperature profiles in the two directions are mostly consistent with each 

other within the error bars and drop to 1 .5 keV at 1 .5 rsoo . As observed by Suzaku in hot clusters, the entropy profile 
tH- ' becomes flatter beyond rsoo, in disagreement with the r'^' relationship that is expected from accretion shock heating 

J> . models. When scaled with the average intracluster medium (ICM) temperature, the entropy profiles of clusters 

^D ' observed with Suzaku are universal and do not depend on system mass. The hydrostatic mass values in the void 

^^ \ and filament directions are in good agreement, and the Navarro, Frenk, and White universal mass profile represents 

the hydrostatic mass distribution up to ^ 2 7'5oo. Beyond rsoo, the ratio of gas mass to hydrostatic mass exceeds 
the result of the Wilkinson microwave anisotropy probe, and at rioo, these ratios in the filament and void directions 
C^ . reach 0.4 and 0.3, respectively. We discuss possible deviations from hydrostatic equilibrium at cluster outskirts. We 

derived radial profiles of the gas-mass-to-light ratio and iron-mass-to-light ratio out to the virial radius. Within rsoo, 
the iron-mass-to-light ratio of the Hydra A cluster was compared with those in other clusters observed with Suzaku. 
Key words: galaxies: clusters: individual (Hydra A cluster) — X-rays: galaxies: clusters — intergalactic 
medium 



.^.' 1. Introduction Imaging Spectrometer (XIS; Koyama et al. 2007), Suzaku 

(Mitsuda et al. 2007) was able to unveil for the first time 

Clusters of galaxies are the largest self-gravitating systems the intracluster medium (ICM) beyond rsoo, a radius within 

in the Universe, and offer unique information on the process of which the mean cluster-mass density is 500 times the cosmic 

structure formation governed by cold dark matter (CDM). In critical density. The accurate calibration of the XIS also al- 

addition, these clusters are considered as a laboratory for study- lows precise measurements of the ICM temperature (Sato et 

ing thermal and chemical evolutions of the Universe in which al. 2011). Suzaku derived the temperature and entropy profiles 

baryons play the most important role. X-ray observations pro- of the ICM of several massive clusters up to the virial radius 

vide valuable information about the structure formation, gas (George et al. 2009; Reiprich et al. 2009; Bautz et al. 2009; 

heating and cooling, and metal enrichment of galaxy clusters. Kawaharada et al. 2010; Hoshino et al. 2010; Simionescu et 

Because the dynamical time-scale of clusters is comparable al. 2011; Akamatsu et al. 2011). From the center to r2oo, a 

to the Hubble time, cluster outskirts should maintain original systematic drop in temperature was found by a factor of ^ 

records of cluster evolution via accretion of gas and substruc- 3, and the observed entropy profiles become flatter beyond 

tures from filaments of the surrounding large-scale structure of rsoo. These profiles are lower than the entropy profile pre- 

the Universe. dieted by the numerical simulations of gravitational collapse 

Thanks to the low and stable background of the X-ray (Tozzi & Norman 2001; Voit, Kay, Bryan 2005), which is pro- 
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portional to r^'^. Kawaharada et al. (2010) discovered that be- 
yond rsoo of the Abell 1689 cluster, the total mass obtained 
from weak-lensing observations with Subaru is larger than that 
calculated assuming hydrostatic equilibrium. Therefore, one 
explanation for the low entropy profiles at cluster outskirts is 
that infalling matter retained some of its kinetic energy in bulk 
motion (Bautz et al. 2009; George et al. 2009; Kawaharada et 
al. 2010). Based on Suzaku observations of the Perseus cluster, 
Simionescu et al. (201 1) proposed a gas-clumping effect as an 
additional interpretation. Hoshino et al. (2010) and Akamatsu 
et al. (2011) discussed possible deviations of electron temper- 
ature from ion temperature to explain the observed lower tem- 
perature and entropy profiles. To clarify these effects in the out- 
skirts, dependence on the system mass or on the ICM temper- 
ature should be examined. With XMM-Newton observations. 
Urban et al. (2011) found a similar flattening of the entropy 
profile in one direction in the Virgo cluster {kT = 2.3 keV). 
However, Humphrey et al. (201 1) detected no evidence of flat 
profiles at large scales (> rsoo) in the relaxed fossil group or in 
the poor cluster, RXJ 1 159H-5531. Therefore, more samples of 
medium-sized clusters are required. 

In addition to measuring the temperature and entropy pro- 
files, Suzaku measured the abundance of Fe in the ICM beyond 
0.5ri8o (Fujita et al. 2008; Tawa et al. 2008; Simionescu et al. 
2011). Metal abundances in the ICM also provide important 
information on the chemical history and evolution of clusters. 
The ASCA satellite first measured the distribution of Fe in the 
ICM (Fukazawa et al. 2000; Finoguenov et al. 2001). Recently, 
spatial distributions of Fe from 0.3 to 0.4 riso have been stud- 
ied with the XMM-Newton and Chandra satellites (Vikhlinin 
et al. 2005; Maughan et al. 2008; Leccardi & Molendi 2008; 
Matsushita 201 1). Since metals are synthesized by supernovae 
(SNe) in galaxies, the ratios of metal mass in the ICM to the to- 
tal light from galaxies in clusters or groups, (i.e., metal-mass- 
to-light ratios) are key parameters in investigating the chemi- 
cal evolution of the ICM. Suzaku measured the iron-mass-to- 
light ratios (IMLR) of several clusters and galaxy groups out 
to 0.2 - 0.5 ri8o (Matsushita 2007; Komiyama et al. 2009; 
Sato et al. 2007; Sato et al. 2008; Sato et al. 2009a; Sato et al. 
2009b; Sato et al. 2010; Sakuma et al. 201 1), and with XMM- 
Newton, the IMLR of the Coma cluster out to 0.5ri8o was de- 
rived (Matsushita et al. 201 1). The IMLR profiles increase with 
radius, indicating that Fe in the ICM extends farther than stars. 

The Hydra A cluster (z = 0.0539) with an ICM tempera- 
ture of ^3 keV is one of the prototype cool-core clusters in 
which McNamara et al (2000) discovered a displacement of 
X-ray gas in the central region through the radio lobes from 
the central active galactic nucleus (AGN). This cluster is also 
known as Abell 780 and remains a major examples of AGN 
interaction that is studied through radio and X-rays (Taylor 
et al. 1990; David et al. 2001; Nulsen et al. 2002; Nulsen et 
al. 2005; Lane et al. 2004; Wise et al. 2007; Simionescu et 
al. 2009a; Simionescu et al. 2009b; Kirkpatrick et al. 2009). 
Through Chandra and XMM-Newton observations, a sharp 
X-ray surface-brightness edge was detected at radii between 
4.3' and 6' (200-300 kpc) and was interpreted as a shock 
wave caused by an AGN outburst. The abundances of the 
Hydra A cluster was measured within ~0.3 riso with XMM 
(Simionescu et al. 2009b; Matsushita 201 1). 
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Fig. 1. XISO image of the Hydra A cluster within 0.5^.0 keV energy 
range. The Ught blue cross indicates the X-ray peak of the cluster. 
Differences in the vignetting effect are uncorrected. Regions of spectral 
accumulation are shown as magenta rings. The ring radii are 2', 4', 6', 
8', 10', 12', 15' and 23'. Black circles represent excluded regions 
around point sources. 

This paper reports the results of two Suzaku observations 
in the northern half of the Hydra A cluster out to 25'(~1.5 
Mpc), which corresponds to the virial radius. The observa- 
tions were conducted using the XIS. This study reports the first 
Suzaku observations of a medium-sized cluster with an aver- 
age temperature of ^3 keV up to the virial radius. The two 
observed fields continue in a filament direction that continues 
to the Abell 754 cluster and a void direction of the large-scale 
structure of the Universe (see Figure 9). 

In this study, we use the Hubble constant, Hq = 
70 km s"^ Mpc^i, f7,„ = 0.27 and Qa = 0.73. The luminos- 
ity distance Z3l and angular size distance Da to the Hydra A 
cluster are I?l — 240 Mpc and Da = 217 Mpc, respectively, 
and 1' corresponds to 66.6 kpc. We use the epoch J2000.0 for 
the right ascension and declination of the equatorial coordi- 
nate system. In addition, we use the solar abundance table by 
Lodders (2003), in which the solar Fe abundance relative to H 
is 2.95 X 10^^. Considering a difference in solar He abundance, 
the Fe abundance yielded by Lodders (2003) is 1 .5 times higher 
than that using the photo-spheric value by Anders & Grevesse 
(1989). Errors were quoted at a 90% confidence level for a 
single parameter 

2. Observation and Data Reduction 

2.1. Suzaku observations 

Suzaku conducted two observations of the Hydra A cluster 
in November 2010, which was during the Suzaku Phase- V pe- 
riod. The details of the observations are summarized in Table 1 . 
The first observation, HYDRA A-1, was 12' north-west offset 
from the X-ray peak of the Hydra A cluster with coordinates 
(RA, Dec) = (139.5236, -12.0955) in degrees. The second ob- 
servation, HYDRA A-2, was 12' north-east offset from the X- 
ray peak of the Hydra A cluster The north-west and north-east 
fields continue into the filament and void structures, respec- 
tively (See figure 9). Hereafter, we refer to the north-west and 
north-east offsets as filament and void, respectively. 
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Table 1. Observation log of the Hydra A cluster 



Field name 



Target name 



Sequence 
number 



Observation 
date 



RA* 

(deg) 



Dec* 

(deg) 



Exposure 

(ks) 



North-west (filament) 
North-east (void) 



HYDRA A- 1 
HYDRA A-2 



805007010 
805008010 



2010-11-08 
2010-11-09 



139.3723 
139.6749 



-11.9472 
-11.9472 



35.6 
34.9 



Coordinates are referred to J2000.0. 



The XIS, which was operated in its normal mode during the 
observations, consists of three sets of X-ray CCD (XISO, XIS 1, 
and XIS3). XISl is a back-illuminated (BI) sensor, while both 
XISO and XIS3 are front-illuminated (FI) sensors. Data re- 
duction was done with HEAsoft version 6.11. The XIS event 
lists created by rev 2.5 pipeline processing were filtered using 
the following additional criteria: a geomagnetic cutoff rigid- 
ity (COR2) >6 GV, and an elevation angle > 10° from the 
earth limb. The 5x5 and 3x3 editing modes data formats were 
added. The exposure times after data selection are shown in 
Table 1. 

To obtain the temperature and electron density profiles, we 
accumulated spectra within annular regions centered on the X- 
ray peak as shown in Figure 1. Regions around calibration 
sources and the innermost region within 2' were excluded from 
the spectral analysis. The reason for this exclusion is that we 
could not obtain sufficient photon statistics to analyze the spec- 
tra (which have multiphase features at the cluster center) be- 
cause the given area was detected on the edge of a CCD chip in 
our observation. In addition, Suzaku's PSF is also insufficient 
to resolve such a complex structure. However, with their pow- 
erful imaging capability, Chandra and XMM observations have 
unveiled such complex features in the central region. Using the 
ewavelet tool in SAS\ we searched for point-like sources in 
the Suzaku images in the energy range of 2.0-10.0 keV. In this 
energy range, the flux levels were approximately > 2 x 10^^^ 
erg cm^^ s^^. These point sources were excluded as circular 
regions with radii of 1' (figure 1). The non X-ray background 
(NXB) was subtracted from each spectrum using a database of 
night Earth observations with the same detector area and COR 
distribution (Tawa et al. 2008). 

We included the degradation of energy resolution due to ra- 
diation damage in the redistribution matrix file (RMF) gener- 
ated by the xisrmfgen Ftools task. In addition, we created 
an ancillary response file (ARE) using the xissimarfgen 
Ftools task (Ishisaki et al. 2007). A decrease in the low-energy 
transmission of the XIS optical blocking filter (ORE) was also 
included in the ARE. For filament and void fields, we generated 
ARE files assuming uniformly extended emission from an en- 
riched region with a 20' radius. We used the XSPEC_vl 2.7.0 
package and ATOMDB_v2.0.1 for spectral analysis. Each 
spectrum was binned and each spectral bin contained a min- 
imum of 50 counts. To avoid systematic uncertainties in the 
background, we ignored energy ranges above 7 keV and below 
0.7 keV. In addition, we excluded the narrow energy band be- 
tween 1.82 and 1.84 keV in the fits because of the incomplete 
response around the Si edge. 
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2.2. XMM observations 

We used the same data and analysis as that reported by 
Matsushita (201 1). The XMM-Newton archival data of Hydra 
A (observation identifier 0109980301) had exposure times of 
14.0 ks, 19.1 ks and 22.2 ks forMOSl, MOS2 and PN, respec- 
tively, after background flares were screened out. We selected 
events with patterns smaller than 5 and 1 3 for the PN and MOS, 
respectively. Spectra were accumulated in concentric annular 
regions of 0'-0.5', 0.5'-l', l'-1.5', 1.5'-2', 2'-3', 3'-4', 4'-5', 
5'-6', 6'-7', 7'-8', and 8'-9', centered on the X-ray peak of the 
Hydra A cluster Here, the X-ray peaks were derived using the 
ewavelet tool of SAS-v8.0.0, and luminous point sources were 
excluded. The spectra from MOSl and MOS2 were added. 
The background spectrum for each annular region was calcu- 
lated by integrating blank-sky data in the same detector region. 
From deep-sky observations with XMM-Newton, we selected 
the data having background most similar to that of the Hydra 
A cluster and the faintest Galactic emission, after screening 
out background flare events from the data and the background, 
following Katayama et al. (2004). Next, we scaled the back- 
ground using a count rate between 10 and 12 keV. 

The response matrix file and the ARE corresponding to each 
spectrum were calculated using SAS-v8.0.0. Further details 
appear in sections 2 and 3 of Matsushita (201 1). 

3. Data Analysis 

3.1. Estimation of background spectra 

To derive the radial profiles of the temperature, electron 
density, and Ee abundance of the ICM, we fitted the NXB- 
subtracted spectra with a thermal plasma model (APEC: Smith 
etal. 2001). 

We first fitted the spectra in the outermost ring (15'-23' re- 
gion) to determine the local X-ray background. As reported 
by Yoshino et al. (2009), the background emission of Suzaku 
XIS can be fitted with a three-component model: two thermal 
plasma models (APEC; Smith et al. 2001) for the local hot 
bubble (LHB) and the solar wind charge exchange (SWCX), 
the Milky Way Halo (MWH), and a power-law model for 
the extragalactic cosmic X-ray background (CXB). MWH and 
CXB components were convolved with an photoelectric ab- 
sorption in the Galaxy, A^h ■ Therefore, we adopted the model 
formula, apecLHB + wabs x (apecMWH + power-lawcxB) as 
background. We assumed a zero redshift and a solar abundance 
for LHB and MWH components. The temperature of the LHB 
and MWH was fixed at 0.1 keV and 0.3 keV, respectively. The 
column density of the Galactic neutral hydrogen was fixed at 
4.7x1020 cm"2 (Kalberla et al. 2005), and the photon index 
of the CXB component was fixed at 1.41 (Kushino et al. 2002). 
Using the formula apecLHB + wabs x {apecicM + apecuwn + 
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Table 2. Resultant parameters of the background components. 



Cosmic X-ray Background Local Hot Bubble MiUcy-Way halo 



normalization* 



normalization^ 



normalization^^ 

xlO-4 



.92: 



-0.85 



.93: 



-5.75 



1.83: 



-1.72 



Reduced-x'^ 

X^/d.o.f 



1.04(144/139) 



* Measured at IkeV with unit of photons cm^'^ s^'^ keV^'^ sr^'^. 

t Normalization of the APEC, component divided by the solid angle, f!'^, assumed in the uniform-sky ARF calculation (20' radius), Norm 
f 71^11 jjdV/{4:TT{l + z)'^D^)/Q" X 10"^'* cm"'' 4007r arcmin^^, where D^j is the angular distance to the source. 



power-lawcxB), we fitted the spectra with a sum of the back- 
ground components and a single-temperature APEC model for 
the ICM. Temperature and normalization of the ICM compo- 
nent were allowed to vary, and redshift was fixed at 0.0539. 
Normalizations of the background components were also left 
free. The spectra of the filament and void fields were fitted 
simultaneously, and each parameter of background and abun- 
dance of the ICM was assumed to be the same. 

Results of the spectral fit are shown in Figure 2. The result- 
ing parameters for the background and ICM are shown in Table 
2 and 3, respectively. The normalizations obtained for the CXB 
and the Galactic components are consistent with those derived 
by Kushino et al. (2002) and Yoshino et al. (2009), respectively. 

Several blank fields observed with Suzaku contained emis- 
sion with kT = 0.6 - 0.8 keV (Yoshino et al. 2009). Therefore, 
we added an APEC component with kT = 0.6 keV and refit- 
ted the outermost spectra. The normalization of the 0.6 keV 
component had a much smaller value than the other Galactic 
components, and the temperatures and normalizations of the 
ICM did not change. 
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Fig. 2. NXB-subtracted spectra of XIS 1 (crosses) for 
the outermost ring (15'-23' region), fitted by the 
apecLHB + wabs X (apecjcM + "-pecMwu + power-lawcys.) 
model (stepped solid lines). Black and red colors con'cspond to the 
directions in the filament and void, respectively. The bottom panel 
shows residuals of the fit. The contributions of the ICM are plotted as 
black and red solid lines, those of the Galactic emission as orange and 
green fines, and those of the CXB as blue lines. 



0) 
Q. 

E 

0) 




filament 

void 
Chandra 

A (keV) b (kpc) 
4.39±0.16 1293±121 
3.72±0.20 1488±339 



50 



100 200 
Radius (l<pc) 



500 1000 



Fig. 4. Radial temperature profiles of the filament (black) and 
void (red) directions. Solid lines show the best fit function of 
A X exp(— r/6). Chandra results (David et al. 2001) are plotted in 
blue. Vertical dotted fines show rsooi^iso. and rioo, derived from 
the fitting of the hydrostatic mass with the NFW universal mass profile 
(Navarro et al. 1997), as described in subsection 4.4. 

3.2. Spectral fits of annular regions 

We fitted the XIS spectra in each region of the Hydra A clus- 
ter using a single-temperature APEC model for the ICM with 
the Galactic absorption, A^h- The spectra of each annular re- 
gion of the two fields were fitted simultaneously by minimiz- 
ing the total x^ value. Here, the common background model 
was included for all regions, where the surface brightness of 
the background components were restricted within the statis- 
tical errors derived for the outermost regions. The tempera- 
ture and normalization of the ICM component were free pa- 
rameters. The abundances in each annular region in the two 
fields were assumed to have the same value. A^h was fixed to 
the Galactic value of 4.7 x lO'^'^ cm^^ in the direction of the 
Hydra A cluster 

The resultant parameters are summarized in Table 3, and the 
best-fit spectra are shown in Figure 3. The spectra are well 
represented with the single-temperature model for the ICM and 
the background. 

4. Results 

4.1. Temperature profile 

Figure 4 shows the radial temperature profiles of the two di- 
rections derived from the spectral fits. At a given radius, the 
temperatures for each direction are consistent within statisti- 
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Fig. 3. NXB-subtracted XISl spectra fitted with the single temperature (APEC) model for the ICM. Black and red colors correspond to directions toward 
the filament and void, respectively. Lower panels in each diagram show residual of fit. The contributions of the ICM are plotted as red and black solid 
Unes, while those of the galactic emission and the CXB are plotted as orange, green and blue lines, respectively. 
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Table 3. Parameters resulting from fitting ICM with the single-temperature APEC model in the energy ranges of 0.7-7.0 keV. 
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Region 



filament 



void 



(arcmin) 



(kpc) 



2' -4' 133 - 266 3.81+^;^ 



kT (keV) abundance (solar) normalization* kT (keV) normalization* 



-0.02 

0^^0+0.01 



-0.14 

3.09+°-2^ 



-0.02 

0.0901H1 
0.035+2:^3 

nni5+°''03 

U.U10_o 0Q2 

fO.0016 

-0.0012 

S+0.0012 

-0.0012 

002^+''°°" 



X /d.o.f. 



4' -6' 
6' -8' 
8' - 10' 
10' - 12' 
12' - 15' 
15' - 23' 



266 - 400 
400 - 533 
533 - 666 
666 - 799 
799 - 997 
997 - 1529 



-0.34 
q ni+0.32 
'J"^^-0.23 
q 11+0.28 
'-'■^^-0.23 

Z.OU_(j 27 

2 '^.A+O-a^ 

^■'J^-0.36 

qf? + 0.33 

1 44+0.47 
^■^^-0.33 



0.33 ±0.10 
0.43 ±0.43 
0.28 ±0.15 
0.18 ±0.13 
0.38 ±0.38 
0.30 ±0.19 
0.10 ±0.08 



0.45-* 



3.18^ 



-0.01 

U.UOO_Q (jog 

U.UZO_Q QQ3 

r,Qio+0.003 

010+°°''^ 

U.UlU_o oQj^ 

0036+''°°^2 



-0.24 
7+0.38 

-0.31 
:;+0.86 

-0.81 

4(5+1.60 

"1+0.33 
-0.31 

1 Qq+0.65 
-'-■•J'^-0.42 



3.07! 
3.161 



1.79! 



0.36-* 



0.0058: 
0.0048^ 



1.22(263/216) 
1.09(146/134) 
1.18(177/150) 
1.18(129/109) 
1.07(99/93) 
1.13(170/150) 
1.19(263/221) 



* Nomiahzation of the APEC, component divided by the solid angle, Q", assumed in the unifomi-sky ARF calculation (20' radius), Norm ■ 
r?ie?i£fdy/(47r(l + z)^D'^)/fl" X 10^^* cm"''' 4007r arcmin"^, where D^ is the angular distance to the source. 



cal errors, except around ^200 kpc from the center The re- 
gion around 200 kpc corresponds to the shock front detected 
by Chandra and XMM in the filament direction (Simionescu et 
al. 2009b; Gitti et al. 2011). From 300 kpc to 1200 kpc, the 
temperatures in each direction decrease with radius in a similar 
manner down to ^ 1 .4 ke V. 

From ROSAT observations, Ikebe et al. (1997) found that the 
spectrum at the central region within 2' exhibits the multiphase 
feature. Although the central 2' region was excluded in our 
spectral analysis, we examined the spectral fits with the two- 
temperature model in regions with r > 2'. The reduced-^^ 
for the two-temperature model was not significantly improved 
in comparison with the single-temperature model. 

To determine the hydrostatic mass of Hydra A, we fitted the 
temperature profile of each direction with the exponential for- 
mula, A X exp(— r/6), where A and b are free parameters with 
units of keV and kpc, respectively, and r is the distance from 
the X-ray peak of the cluster center The resulting fit parame- 
ters and best-fit functions are shown in Figure 4. In the filament 
direction, the normaUzation factor A of the exponential model 
was 18% higher than that in the void direction, which reflects a 
temperature difference of ^ 200 kpc. The parameter b for both 
directions is consistent within statistical error 

4.2. Deprojected electron density profile 

We calculated three-dimensional electron density profiles 
from the normalization of the ICM component derived with 
Suzaku, considering the geometrical volume that contributes 
to each two-dimensional region. Here, we assumed spherical 
symmetry within each direction. The electron density profile 
was also obtained using XMM data. To avoid uncertainties in 
the background, we fitted the annular spectra of XMM MOS 
with the APEC model in an energy range of 1 .6-5.0 keV. Here, 
the temperature and Fe abundance were left free within the er- 
ror bars derived by Matsushita (2011). We fitted the derived 
emissivity profiles with a double-/3-model, and calculated the 
three-dimensional radial profile of the electron-density. 

Figure 5 shows the electron density profile derived from 
XMM and Suzaku (this work), Chandra (David et al. 2001) and 
the ROSAT surface-brightness profile scaled with the XMM 
density profile at 200 kpc (Ikebe et al. 1997). The deprojected 
electron density profile derived from XMM and Suzaku obser- 
vations agree well with those from ROSAT within 500 kpc and 
from Chandra within 200 kpc. 

At a given distance from the X-ray peak, electron density 
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distance (kpc) 

Fig. 5. Deprojected electron density profile of the filament (black) and 
void (red) directions. Black and red solid lines show the best fit with 
the double-/? model. Green solid line shows the electron density profile 
derived with XMM data, the light blue solid line shows that derived 
from the ROSAT surface-brightness profile (Ikebe et al. 1997) scaled 
with the XMM profile, and the blue filled circles show that derived 
from Chandra results (David et al. 2001). 

beyond 400 kpc in the filament direction was systematically 
higher than that in the void direction, reflecting the differ- 
ence in normalizations (Table 3 and Figure 3). At 900 kpc, 
for example, the deprojected electron densities in the fila- 
ment and void directions were (2.54 ± 0.18) x 10~* cm~^ and 
(1.76 ± 0.21) X lO""* cm"^, respectively. The electron density 
dropped to ^lO""' cm~^ around 1200 kpc. 

We fitted the density profiles from XMM and Suzaku results 
with the double-/? model from 40 to 1200 kpc. The results 
are shown in Figure 5, and the resulting parameters are given 
in Table 4. From 400 kpc to 1530 kpc, the radial profiles of 
the electron density are also well represented by a power-law 
profile with slopes of r = 1.60 ±0.05 and 1.72±0.07in the fil- 
ament and void directions, respectively. The slope of the elec- 
tron density in the void direction was slightly steeper than that 
in the filament direction. Fitting the electron density profiles 
beyond 100 kpc with the single /3-model resulted in /3 = 0.61 
± 0.02 and (i = 0.1 A ± 0.03 for the filament and void direc- 
tions, respectively. These values are close to the value of 0.68 
derived from the ROSAT surface-brightness profile (Ikebe et 
al. 1997). 



No. ] Suzaku Observations of the Hydra A Cluster 

Table 4. Parameters resulting from fitting tlie deprojected density profile derived from XMM and Suzaku with the double-/3 model*. 



Field 



narrower component 



wider component 



710,1 (10 cm ■=') rc,i (kpc) 



/3i 



no,2 (10 -^cm ■^) rc.2 (kpc) 



/32 



void 
filament 



6.73 ±0.17 
6.61 ±0.19 



21.3±4.7 
22.3±5.4 



0.59 ±0.01 
0.62 ±0.01 



2.16±0.14 
2.73±0.15 



278 ± 264 
193 ±114 



0.95 ±0.07 
0.60 ±0.03 



* double-/3 model is given by the form n{r) = no,i(l + (r/rc^i)'^) '^'^i/'^ + no, 2(1 + (''"/f'c,2)'^) 
The errors of no,i and /3i are calculated by fixing the rc,i and wider component, and vice versa. 



-3^/2 



Table 5. Resultant fit pai'ameters of the deprojected density profile derived from XMM and Suzaku with the single-/? model* beyond 100 kpc. 



Field no (lO"'' cm"^) 



(kpc) 



(3 



void 7.38 ±0.42 160 ±13 0.74 ±0.03 

filament 7.98 ±0.50 131 ± 10 0.61 ±0.02 

' single-/? model is given by the form n{r) = no(l + {r j r c)^')^'^^ ^ '^ . 
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Fig. 6. Entropy profiles in the filament (black) and void (red) direc- 
tions. Crosses and diamonds show entropy calculated from the result- 
ing temperature and electron density, respectively, as shown in Figures 
4 and 5. Chandra results (David et al. 2001) are represented by blue 
crosses. The green line shows the entropy profile calculated by combin- 
ing the electron density from XMM and the temperature from Suzaku. 
The light blue line extrapolates the Suzaku and XMM data (green line) 
with a power-law formula fit for the data beyond 200 kpc, which has a 
fixed index of 1 . 1 . 

4. 3. Entropy profiles 

Entropy, which is a useful tool for investigating the thermo- 
dynamic history of hot gas in groups and clusters, is defined 
as 

kT 



S^ 



2/3' 



(1) 



where T and rie are the temperature and deprojected electron 
density, respectively. In Figure 6, we show the entropy pro- 
files of the Hydra A cluster calculated from the derived tem- 
perature and electron density profiles shown in Figures 4 and 
5. Entropy in the two directions increased from the center to 
800 kpc, or rsoo, which was derived from the Navarro, Frenk 
and White(NFW) model fit to the hydrostatic mass described in 
section 4.4. Between 500 kpc and 800 kpc, the entropy in the 
void direction was systematically higher than in the filament 



direction. At the two outermost radial bins (beyond rsoo), the 
entropies of the two directions agree well. In addition, we cal- 
culated the entropy profile by combining the parameters result- 
ing from the fit of the electron density profile from XMM and 
the temperature profile from Suzaku as indicated by the green 
line in Figure 6. The entropy profiles using Suzaku, XMM, and 
Chandra data are consistent with each other (figure 6). 

The derived entropy profiles were compared with a power- 
law model with a fixed index of 1.1, which was expected from 
simulations of the accretion shock heating model (Tozzi & 
Norman 2001; Ponman, Sanderson, Finoguenov 2003). From 
200 kpc to 500 kpc, the derived entropy profiles agreed well 
with the r^^ relationship. In contrast, beyond 800 kpc, or 
^ rsoo, the entropy profiles in the two directions were system- 
atically lower than the predicted by the r^^ relationship. Inside 
200 kpc, the entropy exceeded the model prediction, which at- 
tribute to the presence of the cool core boundary (Peres et al. 
1998) and the shock front around 200 kpc. 

4.4. Hydrostatic mass 

We estimated the hydrostatic or gravitational mass assuming 
the hydrostatic equilibrium given by, 



MH.E.(<r) 



kT{r)r /dlnp{r) dlnkT{r) 



dlnr 



din? 



(2) 



where A/h.e.(< r) is the hydrostatic mass within the three- 
dimensional radius r, G is the gravitational constant, k is the 
Boltzmann constant, p is the gas density, kT is the tempera- 
ture, and yUJTip is the mean particle mass of the gas (the mean 
molecular weight is /i = 0.62). Here, we assumed spherically 
symmetric mass distribution. 

On the basis of the temperature and electron density pro- 
files, we calculated the hydrostatic mass of the Hydra A clus- 
ter beyond 100 kpc. The temperature gradient of the ith shell 
was calculated from a power-law fit of temperatures of shells 
(i — 1), i, and [i + 1). We employed the electron density pro- 
file from the single-/? model fit shown in table 5. The derived 
hydrostatic mass is plotted in Figure 7 as diamonds with the 
integrated gas mass profiles and stellar mass profiles derived in 
subsection 4.6. The derived mass was calculated for the whole 
azimuthal angle under the assumption of spherical symmetry 
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Fig. 7. left:Integrated hydrostatic and gas-mass profiles. All the black and red colors show the filament and void directions, respectively. The masses 
were calculated in the whole azimuthal angle under the assumption of spherical symmetry for each direction. The hydrostatic mass was calculated from 
the temperatures of shells (« — 1), i, and (i + 1) (diamonds) and the exponential relationship (dashed lines). All dashed lines show 90% upper and lower 
Umits. The green line shows the best-fit NFW model. The integrated stellar mass profiles in the filament and void directions using 2MASS data are 
represented by stars. Blue stars are the integrated stellar mass profile in the whole azimuthal angle coiTected for the contribution of the fainter galaxies 
using LF (solid). right:Cumulative profiles of gas fraction, ,fgas(<'') assuming hydrostatic equilibrium in both directions. The diamonds show /gas(< ^) 
of the integrated hydrostatic and gas mass of the corresponding directions. Dashed lines represent /gas (< r) applying the result of the NFW model fit. 
The value of the cosmic baryon fraction (light blue) from WMAP (Komatsu et al. 2011) is plotted. The stars represent the stellar fraction, /star(< t), 
using the NFW model fit. The meanings of the colors and fine styles are the same as in the left panel. 



using the parameters derived for each direction: thus, the mass 
in each direction was four times larger than that in the corre- 
sponding direction. Although the electron density in the fil- 
ament direction at a given radius was higher than that in the 
void direction, the difference in the slopes of the electron den- 
sity was small. As a result, the hydrostatic masses derived for 
the two directions agree well. In addition, we calculated the 
90% range of hydrostatic mass using the exponential relation- 
ship of the temperature profiles and fits of the electron-density 
profiles with the double-/? model, which are plotted as dashed 
lines in Figure 7. 

We fitted the derived hydrostatic mass profiles with the fol- 
lowing NFW formula, which shows the equilibrium density 
profile of the dark matter halo (Navarro et al. 1996; Navarro 
et al. 1997); 



5c = 



200 



(4) 



MNFw(<r)=47r(5cPcrs 



ln(l + a;) 



1 



,x = r/rs,(3) 



where AfNFw(< ^) is the mass within a radius r, p^ is the 
critical baryon density of the Universe, r^ is the scaled radius, 
and (5c is the characteristic density, that can be expressed in 
terms of the concentration parameter (c = tiqqJTs) as 



3 ln(l + c)-c/l + c' 

Since the hydrostatic masses derived for the two directions 
agree well, we fitted the two hydrostatic mass profiles with 
Suzaku and that with Chandra simultaneously, excluding the 
shock region from 100 kpc to 200 kpc. The parameters derived 
from the fit are shown in Table 7. The quantity rioo is equal 
to Ttjir in Hydra A cluster, since for our fiducial cosmologi- 
cal model, A„ir — 100 for halos with a redshift z ~ 0.0539 
(Nakamura & Suto 1997). We estimated r^ for various over- 
densities: A = 500, 200, 180 and 100, and ta derived from 
fits with the NFW model are 811 ± 10 kpc, 1189 ± 96 kpc, 
1243 ± 16 kpc, and 1577 ± 22 kpc, respectively. The rioo is 
close to the outermost ring radius. Possible deviations from hy- 
drostatic equilibrium have been previously discussed to explain 
the entropy flattening beyond r5oo (Bautz et al. 2009; George 
et al. 2009; Kawaharada et al. 2010): therefore, we fitted the 
hydrostatic mass out to tsoq. We were able to obtain almost 
the same results as shown in Table 7. Thus, the NFW model 
represents the hydrostatic mass out to rioo within error bars as 
shown in Figure 7. 
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Fig. 8. H.E. mass density profiles for ttie filament (black) and void 
(red) directions. Density was calculated from the integrated H.E. mass 
profile. The r~^ and r~^ profiles are plotted with solid blue and or- 
ange lines, respectively. 

Moreover, we calculated the total mass density profile from 
the total integrated hydrostatic mass profile in each direction. 
The density profile was also calculated from the NFW model 
obtained above. The calculated density profiles were compared 
with r^^ and r^^ profiles (Figure 8). The NFW density profile 
varied from r^^ to r^^ with increasing radius. Within r^Qo, 
the densities are well represented by the r^^ profile, which is 
shallower than the asymptotic matter density slope p ex r~^. 
Beyond rsoo, the r^^ profile is preferable than the r^^ profile. 

4.5. Gas fraction 

We derived the cumulative gas mass fraction as. 



/gas(<r) 



Mgas(<r) 

A/H.E.(<r)^ 



(5) 



where A/gas(< r) and A/h.e(< r) are the gas mass and hy- 
drostatic mass, respectively, within a sphere of radius r. As 
shown in Figure 7 and Table 6, the gas fraction increased in 
the outskirts and reached the cosmic mean baryon fraction de- 
rived from seven-year data of Wilkinson microwave anisotropy 
probe (WMAP7; Komatsu et al. 2011) at -- 0.5 tsoq. At 
'^looi /gas(< ?') in the filament direction exceeded the result 
of WMAP7, while the lower limit of /gas(< r) in the void 
direction reached the WMAP 7 fraction. Adopting the fit of 
the NFW model to the hydrostatic mass resulted in /gas(< r) 
in the void and the filament directions at rioo, exceeding the 
cosmic baryon fraction by a factor of 2 and 3, respectively. 
Furthermore, /gas(< r) in the void direction at rigg was close 
to that derived for the north-west direction of the Perseus clus- 
ter at r2oo (Simionescu et al. 2011): rigo was derived from the 
NFW model fit of the hydrostatic mass. The radial profile of 
/gas(< r) of the Abell 1246 cluster with kT = 6.0 keV exhib- 
ited similar behavior as that of the our result of the Hydra A 
cluster (Sato et al. 2012). 

4.6. K-band Luminosity of galaxies and stellar fraction 

The ratio of gas mass to stellar luminosity is a key pa- 
rameter for studying the star formation efficiency in clusters 



of galaxies. Since almost all metals in the ICM are synthe- 
sized in galaxies, the metal-mass-to-hght ratio provides a use- 
ful measure for studying the ICM chemical evolution. To es- 
timate the gas-mass and Fe-mass-to-light ratios, we calculated 
the A' -band luminosity of Hydra A cluster on the basis of the 
Two Micron All Sky Survey (2MASS) catalog, whereby all 
data is presented in a 2°x2° box: the coordinates of objects 
are shown in the left panel of Figure 9. The smoothed dis- 
tribution of galaxies around Hydra A obtained from 2MASS 
is shown in the right panel in Figure 9. The distribution of 
galaxies is elongated in the filament direction toward the Abell 
754 cluster, which exhibited a redshift close to that of Hydra 
A. In the void direction, the density of galaxies was signif- 
icantly lower. Toward the south, an additional filament-like 
distribution of galaxies appeared from the Hydra A cluster 
Hydra A itself was discovered to have an apparent magnitude 
of rriK = 11.07, or Lk = 5 x lO^^iif^© using a luminosity dis- 
tance Dl ~ 240 Mpc and a foreground Galactic extinction of 
Ak = 0.015 (Schlegel, Finkbeiner, Davis 1998). Among galax- 
ies in the NED data base^ with available redshifts, those with 
heliocentric velocities differ from that of the Hydra A cluster 
by >2000 km/s were excluded from the analysis. The /v -band 
surface brightness profile centered on Hydra A is shown in fig- 
ure 10. Beyond ^2000 kpc from the cluster center, the AT-band 
surface brightness flattened. Therefore, we adopted the region 
between 2000 kpc and 4000 kpc as a background. 

Figure 10 also shows the luminosity function (LF) of galax- 
ies in the Hydra A cluster within rioo- Here, we used as back- 
ground the LF of galaxies from 2000 kpc to 4000 kpc from the 
center of the Hydra A cluster The completeness limit of galax- 
ies of 2MASS is Kg = 13.5 in apparent magnitude ^, which 
corresponds to the i^-band luminosity, Lk=7 x lO^^L^.o for 
the distance of the Hydra A cluster The total A'-band lumi- 
nosity of galaxies above the 2MASS complete limit and that 
of all the detected galaxies including those below the limit are 
1.5 xIO^^Lk^o and 2.5 x 10^^ A^.o, respectively. To calcu- 
late the total light of the cluster, we evaluate the contribution 
of fainter galaxies below the 2MASS limit. The LF of galaxies 
in clusters can be described with the Schechter function as. 



4){L)dL = (f)^ 



t) -p(-vi.)<^ 



(6) 



where the parameters L», 0,^, and a are the characteristic lu- 
minosity, number densities, and faint-end power-law index, re- 
spectively. Upon adopting the A* of stacked cluster galaxies 
by Lin, Mohr, Stanford (2003) and integrating the Schechter 
function below the 2MASS limit, the total A'-band luminos- 
ity becomes 2.2 xlO^^A^.o, which is smaller than that of 
all the galaxies detected by 2MASS including those below 
the complete limit. If we allow A, to vary, the fit for the 
LF of the Hydra A cluster improves, as shown in Figure 10, 
and the total K-band luminosity becomes 2.8 xlO^^Ai^-©. 
Therefore, the simple sum of AT-band luminosities of the de- 
tected galaxies, including those below the 2MASS complete 
limit is close to that integrated luminosity of the Schechter 
function. Therefore, we adopted the total AT -band luminosity 



^ http://nedwww.ipac.caltech.edu/ 

^ See http://www.ipac.caltech.edu/2mass/releases/second/doc/explsup.html 
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Table 6. Cumulative values of the hydrostatic mass (Mh.e.), gas mass (Mgaa), and gas fraction (fgas) of the Hydra Hydra A cluster. 



[Vol. , 



Region 




filament 






void 




r 
(kpc) 


Mgas 

(10'-* Mq) 


Mh.e. 
(10^* Mq) 


Jgas 


Mgaa 
(10" Mq) 


Mh.e. 

(10" Mq) 


J gas 


133 - 266 


0.052 ±0.001 


0.39 ±0.04 


0.13 ±0.01 


0.052 ±0.001 


0.36 ±0.02 


0.15 ±0.01 


266 - 400 


0.11 ±0.01 


0.72 ±0.07 


0.15 ±0.02 


0.11 ±0.01 


0.68 ±0.06 


0.16 ±0.02 


400 - 533 
533 - 666 


0.19 ±0.02 
0.26 ±0.03 


1.06 ±0.09 
1.27 ±0.21 


0.17±0.02 
0.21 ±0.04 


0.17 ±0.02 
0.24 ±0.03 


1.03 ±0.13 
1.43 ±0.39 


0.17 ±0.03 
0.17 ±0.05 


666 - 799 
799 - 997 


0.35 ±0.05 
0.46 ±0.07 


1.35 ±0.35 
2.04 ±0.37 


0.26 ±0.08 
0.22 ±0.05 


0.30 ±0.04 
0.38 ±0.06 


3.03 ±1.41 
1.66 ±0.31 


0.10 ±0.05 
0.23 ±0.06 


997-1529 


0.71 ±0.12 


1.82 ±0.60 


0.39 ±0.14 


0.55 ±0.11 


1.67 ±0.78 


0.33 ±0.17 



Table 7. Resulting radii Tc and r2oo. and concentration parameter C200 derived from the NFW model. 



(kpc) r2oo (kpc) 



fitted up to rsoo 
fitted up to rioo 



126 ±8 
123 ±8 



1189 ±! 
1183 ±< 



9.44 ±0.48 
9.62 ±0.49 



of (2.5 ±0.3) xlQi^i^Q. 

We collected A' -band luminosity of galaxies in the north- 
west and north-east sectors for the filament and void direc- 
tions, respectively, as well as in the whole azimuthal angle and 
derived brightness profiles in the /f-band. For the whole az- 
imuthal angle, we collected galaxies above the 2MASS com- 
plete limit and evaluated the contribution of fainter galaxies 
by integrating the Schechter function. We obtained almost 
the same radial profile of A'-band luminosity as that from the 
simple sum of luminosities of all the galaxies detected with 
2MASS. Therefore, for the filament and void directions, we 
collected all the galaxies detected by 2MASS in each direction 
and derived radial profile of K-band luminosity. After subtract- 
ing the background, we deprojected the brightness profile of 
galaxies (under the assumption of spherical symmetry) and de- 
rived three-dimensional radial profiles of A'-band luminosity 
in the filament, void, and whole azimuthal angle. 

Figure 7 shows the integrated AT-band luminosity profiles. 
As opposed to the integrated luminosity of the whole azimuthal 
angle, figure 7 uses the four-times luminosity of the filament 
and void directions. The filament direction shows a system- 
atically higher A'-band luminosity than that in void direction. 
The integrated A' -band luminosity profiles of these two direc- 
tions become flatter from ~^800-1000 kpc, and the contribution 
of galaxies beyond 1000 kpc is small. The integrated luminos- 
ity profile of the whole azimuthal angle is between those in the 
filament and void directions except for the outermost region 
around 1200 kpc (Figure 7 and 9). This result is attributed to 
the luminosity of the whole azimuthal angle being affected by 
the filamentary structure around 1200 kpc in the south of the 
Hydra-A cluster 

We derived cumulative stellar mass fraction as. 



/star(< r) 



Mstar(< r) 



(7) 



MH.E.(<r) 

where, Mstar(< f) is the total stellar mass within r, derived 
from the total K-band luminosity within r, under the assump- 
tion that the stellar-mass-to-light ratio in the K-band is unity 
(Nagino & Matsushita 2009). At a given radius, the /star(< r) 



in the filament direction is higher by a factor of 2-3 than that 
in the void direction, reflecting the concentration of galaxies in 
the filament direction (Figure 7). Beyond r5oo, ./star(< r) is 
smaller than ^ 1%, and does not contribute to the baryon frac- 
tion. We also derived the stellar mass fraction using the whole 
azimuthal angle using the NFW mass for the filament and void 
directions. Although the stellar fraction continues to decrease 
with radius in the filament and the void directions, the stel- 
lar fraction in whole azimuthal directions shows a minimum at 
~ ^500, and increases with radius by a factor of '^ 1 .5 from rsoo 
tor 100. 

Using the integrated AT-band luminosity profile, we calcu- 
lated the gas-mass-to-light ratio in the filament and void direc- 
tions (Figure 11). The radial profile of the gas-mass-to-light 
ratio in the two directions increased with radius out to rioo. 
We also calculated the gas-mass-to-light ratio in the whole az- 
imuthal angle, using the stellar luminosity in the whole az- 
imuthal angle and the gas mass in the void direction which 
agrees that for the whole azimuthal angle derived by ROSAT 
(Ikebe et al. 1997). 

5. Discussions 

The Hydra A cluster is the first medium-sized cluster ob- 
served with Suzaku out to the virial radius. Suzaku observed 
the northern half of the cluster, and we derived radial profiles 
of temperature, electron density, entropy, hydrostatic mass, and 
K-band stellar luminosity in the filament and void directions. 
Temperature and entropy carry important information about 
the thermal history of the ICM. In subsection 5.1, we com- 
pare these profiles of the Hydra A cluster and other clusters 
observed with Suzaku and XMM to study their dependence 
on the system mass or on average ICM temperature. Next, 
we discuss possible explanations for entropy flattening and the 
higher gas-mass-to-hydrostatic-mass ratio at cluster outskirts. 
Discrepancies between electron and ion temperatures are dis- 
cussed in subsection 5.2, deviations from hydrostatic equilib- 
rium and the gas-clumping effect are discussed in subsections 
5.3 and 5.4. We also derive IMLR of the Hydra A cluster out 



No.] 



Suzaku Observations of the Hydra A Cluster 



11 



o 
'■*—' 

CO 

"o 
Q 





. \ ' ' ' ' \ 




1 1 1 J 1 1 


' Jf' 


' 


y- 






/ 


- • 


m 


• 






m 




: •* 






■ 


11.5 


: t. i^ 




• 


•. 


•• 


• 

* * - 




. • 




• •. 


• . 


• 


• - 






• 


-^H-^ 


• 


d* 


• 
• •• 


-12 


■ • . 










• 




•■■.* 




'^S 


• • 

• 
• 

• 




• 


12.5 


• 




. :, * 




_ 








•*•: 


• 








• 




*•/* 






• 
• • 










• 




• 


-13 


* • 




, , , -1 , , 


. . 1 


• 


• 
. . 1 



Hydra-A ( 25.0 [sqdeg] ) 



. »^ 



140.5 140 139.5 139 

Right Ascention (deg) 



138.5 



Ua 141.5 141 140.5 140 139.5 139 130.5 138 137.5 
BA [deg] 



Fig. 9. Two-dimensional density distributions of galaxies around Hydra A, retrieved from 2MASS in the ftT-band. left: Stars represent galaxies with 
available redshift in the NASA/IPAC Extragalactic Database (NED). Green and blue stars indicate members and non-members, respectively, of the Hydra 
A cluster based on their redshift. The boxes represent the Suzaku fields of view (FOVs) of XIS pointings. The directions of the filament (north-west) and 
the void (north-east) are separated as dashed lines. Right: The box size is 5° x5°. Abell 754 cluster, with coordinates (RA, Dec) = (137.2087, -9.6366) 
in degrees, is shown in top right. A Gaussian smoothing scale is represented by the white hatched circle at bottom right. 
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Fig. 10. left:Surface-brightness profile of Hydra A observed with 2MASS photometric data in the A'-band. Blue points denote the detected galaxies of 
Hydra A and red stars denote the background-subtracted data. The dashed Une shows the subtracted background level and the dotted line shows rioo. 
right: The LF of Hydra A observed with 2MASS photoinetric data in the A'-band within rioo in absolute magnitude form, including background galaxies 
(black diamonds) and that the subtracted background contribution (red crosses). The vertical dashed line represents the complete limit of 2MASS. Blue 
dashed Une and green solid hne show the best-fit Schechter functions when L* was fixed to the result from the stacked LF of galaxies in clusters by Lin, 
Mohr, Stanford (2003) and allowed to vary, respectively. 



to the virial radius in subsection 5.5 to constrain the star for- 
mation history in the cluster 

5.1. Comparison of temperature and entropy profiles with 
other systems 

Figure 12 compares the temperature profiles of the Hydra 
A cluster and several clusters observed with Suzaku, scaled 
with average temperature, < kT >, and ?'2oo- Here, r2oo = 
2.77/iTO yJ{kT)ll{) kcV /E{z) Mpc, with E{z) = (17^(1 + 
z)3 + 1 - r2ji/)i/2 (Henry et al. 2009). Then, raoo of the Hydra 
A cluster is 1.48 Mpc. Within 0.5r2oo, the temperature profiles 
agree with the scaled temperature relationship by Pratt et al. 



(2007), which is derived from a sample of 15 clusters observed 
with XMM out to ^ r^oQ. In contrast, beyond 0.57'200i the ob- 
served temperatures scatter by a factor of two, and tend to be 
lower than the XMM relationship. In the scaled temperature 
profiles, there was no evidence of any significant dependence 
on the average ICM temperature, or system mass. Some vari- 
ation in the temperature profiles should reflect the azimuthal 
variation as observed in Abell 1689 (Kawaharada et al. 2010) 
and Abell 1246 (Sato et al. 2012). 

In figure 13, we compare the entropy profile of the Hydra 
A cluster with other clusters observed with Suzaku; these pro- 
files are scaled with the average ICM temperature and r2oo by 
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Fig. 11. Cumulative gas-mass-to-light ratio Mgas/ Lk in the filament 
(black) and void (red) directions and for the whole azimuthal angle 
(blue). 
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Fig. 12. Radial profiles of temperature of clusters observed with 
Suzaku and XMM. Here, the temperature profiles are scaled with 
the average ICM temperature, < kT >, and are azimuthally aver- 
aged. We plotted the radial profiles of cluster of galaxies Abell 1689 
(Kawaharada et al. 2010), Abell 2142 (Akamatsu et al. 2011), Abell 
1413 (Hoshino et al. 2010), PKS-0745 (George et al. 2009), the Perseus 
cluster (Simionescu et al. 2011), Abell 1246 (Sato et al. 2012), the 
Virgo cluster (Urban et al. 2011), RXJ1159 (Humphrey et al. 2011) 
and the Hydra A cluster (current study). The scaling radius r2oo was 
calculated by Henry et al. (2009). The solid line represents the scaled 
temperature profile by Pratt et al. (2007). 

Henry et al. (2009). Here, we used the weighted average of the 
radial profiles of different directions of each cluster Contrary 
to the expected r^-^ relationship, these profiles become flat be- 
yond 0.5r2oo. Except for PKS-0745 (George et al. 2009), the 
scaled entropy profiles of these clusters were universal with 
ICM temperatures above 3 keV and did not depend on the ICM 
temperature. For the PKS-0745 cluster, Eckert et al. (2011) 
determined that the ROSAT surface-brightness profile is sta- 
tistically inconsistent (7.7cr) with Suzaku results beyond r2oo 
and that the difference is likely explained by the existence of 
additional foreground components at the low Galactic latitude 
of the source. These components were not considered in the 
Suzaku background modeling. 
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Fig. 13. Radial profiles of entropy of cluster observed with Suzaku 
and scaled: with the average ICM temperature, < kT >, and are az- 
imuthally averaged. Data for Abell 1689 (Kawaharada et al. 2010) are 
plotted in dark blue, for Abell 2142 (Akamatsu et al. 2011) in green, 
for Abell 1413 (Hoshino et al. 2010) in orange, for PKS-0745 (George 
et al. 2009) in puiple, for Perseus (Simionescu et al. 201 1) in light blue, 
for Abell 1246 (Sato et al. 2012) in red, and for Hydra A (current study) 
in black. The scaling radius r200 was calculated by Henry et al. (2009). 
The solid line shows the power-law profile with slope of 1 . 1 . The nor- 
malization of the power-law is scaled with our result for the Hydra A 
cluster. 

5.2. Ion-electron relaxation in the HydraA cluster 

One interpretation for the flattening of the entropy profile 
is that caused by deviations in electron and ion temperatures 
(Hoshino et al. 2010; Akamatsu et al. 2011). The time-scale 
for thermal equilibration between electrons and ions through 
Coulomb scattering is given by 

/ n \^^ f kT \ "^^^ 



0.14 



(Takizawa & Mineshige 1998; Takizawa 1998; Takizawa 1999; 
Akahori & Yoshikawa 2010). Considering this time-scale and 
the shock propagation speed, the radial length above which the 
electron temperature is significantly lower than the ion tem- 
perature is proportional to the square of the ICM temperature 
(Takizawa 1999). On the other hand, r2oo is proportional to the 
square root of the ICM temperature. Therefore, some depen- 
dence on the mean ICM temperature is expected in the tem- 
perature and entropy profiles. However, there is no system- 
atic dependence on the average ICM temperature in the tem- 
perature and entropy profiles, as shown in figure 12 and figure 
13. Furthermore, if the ICM temperature is underestimated and 
the actual entropy profile of the Hydra A cluster follows the 
7'!^ relationship, a flat ICM temperature profile is necessary, 
i.e., r~" out to rioo, considering the electron density profile 
of ^-(i 6~i.7) beyond 400 kpc. However, the flat temperature 
profile contradicts to the results of numerical hydrodynamical 
simulations, where negative gradients of the temperature pro- 
files outside cool core regions are naturally produced (Borgani 
et al. 2002; Burns et al. 2010; Nagai & Lau 2011). Therefore, 
for the Hydra A cluster, it is difficult to explain the flattening 
of the observed entropy profiles by deviations of ion-electron 
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temperatures. 

5.3. Deviations from hydrostatic equilibrium 

The remaining interpretation for the flattening of the en- 
tropy profiles at the cluster outskirt is deviations from hydro- 
static equilibrium, because when the gas-clumping effect is 
significant, we also expect that the ICM also deviates from 
hydrostatic equilibrium. Then, infalling matter may have re- 
tained some of its kinetic energy in the form of bulk motions, 
and the thermal energy deficit in the ICM yields lower en- 
tropy. Furthermore, the underestimate of the gravitational mass 
and/or gas clumping leads to an overestimation of the baryon 
fraction. Recent numerical simulations by Vazza et al. (2009) 
and Nagai, Vikhlinin, Kravtsov (2007) have shown that the ki- 
netic energy of bulk motion carries ^ 30% of the total energy 
around the virial radius. With simulations, Nagai & Lau (201 1) 
showed that beyond r2oo, gas clumping leads to an overestima- 
tion of the observed gas density and causes flattening of the 
entropy profile. If the hydrostatic mass was underestimated by 
130% in the filament direction and 100% in the void direction, 
the baryon fractions at 7'ioo would have the same value as that 
of WMAP7, although the observed gas fraction calculated from 
the hydrostatic mass in the void direction is consistent with the 
WMAP7 result within error bars. 

The distributions of the ICM and galaxies in the Hydra A 
cluster elongate in the filament direction. If the ICM is un- 
der hydrostatic equilibrium in a non-spherical dark matter halo, 
the X-ray emission and dark matter halo elongate in the same 
direction. Numerical simulations by Burns et al. (2010) de- 
termined that gas, galaxies, and dark matter continue in the 
filament direction due to the accretion of gas and subclusters. 
Then, it is reasonable that the dark matter halo also elongates 
in the filament direction. However, the hydrostatic mass in the 
void and filament directions of the Hydra A cluster are in good 
agreement. These results indicate that deviations from hydro- 
static equilibrium should be more significant in the filament 
direction than in the void direction. 

The clumping effect will also be higher in the filament di- 
rection, since smaller systems are thought to be accreted from 
this direction. For the Hydra A cluster, we were not able to 
detect clump candidates with Suzaku. The spectra of detected 
point sources are consistent with those of background active 
galactic nuclei. Considering the luminosities of detected point 
sources, the luminosity of clumps should be smaller than a few 
lO''^ erg s^^, which corresponds to cores of groups of galax- 
ies. Therefore, if clumps are significant, their scales should be 
smaller than those of small groups of galaxies. 

5.4. The concentration-mass relationship 

From numerical simulations, a weak variation in the concen- 
tration is expected from low-mass to high-mass clusters, re- 
flecting differences in the formation epochs of low-mass and 
high-mass halos (Navarro et al. 1997; Bullock et al. 2001; 
Dolag et al. 2004). Pointecouteau, Amaud, Pratt (2005), Ettori 
et al. (2011) and Okabe et al. (2010) investigated the relation- 
ship between the concentration parameter and the cluster mass, 
or the concentration-mass (c — M) relationship. The results 
are compatible with the intrinsic dispersion of theoretical pre- 
dictions. Figure 14 compares the c — M relationship of the 
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Fig. 14. Observed distribution of concentration parameters as a func- 
tion of cluster masses. Tlie mass of weak lensing (blue) and X-ray (light 
blue) were defined by Okabe et al. (2010) and Ettori et al. (2010);Ettori 
et al. (2011), respectively. Okabe et al. (2010) used A^no and chq, 
while Ettori et al. (2010);Ettori et al. (2011) used M200 and C2oo. We 
considered J\/ioo for Hydra A. The solid line indicates the median re- 
lationship determined through CDM simulations for the WMAP5 cos- 
mological model, while the region enclosed within the dashed fines cor- 
responds to tfie range of o-(logioc) = 0.1, witfiin wfiicfi most of sim- 
ulated clusters are distributed (Duffy et al. 2008). The large diamond 
shows the observed mass profile of Hydra A (cun'ent work) assuming 
hydrostatic equilibrium. Black and red open triangles correspond to 
cases in which the /gas(< '") value of Hydra A (current work) is the 
same as the cosmic mean baryon fraction (Komatsu et al. 20 11) in the 
filament and the void directions, respectively. 

Hydra A cluster and other clusters observed with weak-lensing 
(Okabe et al. 2010), and clusters observed with XMM by Ettori 
et al. (201 1). An assumption of cluster dynamical state or hy- 
drostatic equilibrium is not required in the weak-lensing re- 
sults. The distribution of the c — M relationship of these clus- 
ters agrees well with the relation expected from the numerical 
simulations by Duffy et al. (2008), although Ettori et al. (201 1) 
found that relatively low-mass systems tend to have higher con- 
centration parameters. 

The hydrostatic mass in the Hydra A cluster of the filament 
and void directions are represented with the same NFW mass 
model with c^i, = 9-10 and Mvir = (2.2 ± 0.4) x 10^^ Mq. 
As shown in Figure 14, this Cvir — A/vii relationship of the 
Hydra A cluster was outside cr(log]^Q c) — 0.1 within which 
most clusters from numerical simulations are distributed Duffy 
et al. (2008). If the hydrostatic mass was underestimated and 
the baryon fraction using the gravitational mass beyond rsoo 
of the Hydra A cluster is the same as the cosmic mean baryon 
fraction in the two observed directions, the fits with the NFW 
model yielded smaller values of Cvii and higher values A/vir as 
shown in Figure 14. In this case, the derived value of A/vir in 
the filament direction became higher than that in the void direc- 
tion. This result also indicates that the ICM in the outskirts of 
Hydra A deviates from hydrostatic equilibrium, and the X-ray 
emission and dark matter halo elongate in the same direction. 
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Fig. 15. Compaiison of iron-mass-to-light-ratio for other clusters. The 
radius is scaled with the characteristic radius r2()0 (Henry et al. 2009). 
For comparison, we plotted the radial profile of cluster of galaxies 
AWM 7 (Sato et al. 2008, Sato et al. in preparation), Abell 262 
(Sato et al. 2009b), the Centaurus cluster (Sakuma et al. 201 1) and the 
Coma cluster (Matsushita et al. 2011) and galaxy groups NGC 5044 
(Komiyama et al. 2009) and NGC 1550 (Sato et al. 2010). Open cir- 
cles are the result from XMM-Newton satellite, and others are Suzaku 
satellite. 

5.5. Metal-Mass-to-Light Ratios out to the virial radius 

We derived integrated IMLR, mass-to-light ratio for Fe, out 
to the virial radius for the first time. Here, we show the IMLR 
profile for the whole azimuthal angle to compare other clusters 
observed with Suzaku. Within 0.3 r2oo and beyond 0.3 r2oo, 
we used the Fe abundances obtained by Matsushita (20 11) and 
this work, respectively. Here, within 8'-15', or 0.35-0.65 r2oo, 
we used the weighted average of the Fe abundance at the three 
radial bins in this region. The error bars of the mass-to-light 
ratio include the errors in both the abundance and the K-band 
luminosity caused by the uncertainty in the contribution of 
fainter galaxies below the 2MASS threshold (see section 5.5 
for details). Figure 15 shows the cumulative IMLR profile of 
the Hydra A cluster The IMLR profile increases with radius 
within 0.4-0.5 r2oo, and becomes flatter beyond 0.5 r2oo with 
the value of ^ 0.01 Mq/ Lk,q, which reflects the flattening of 
the gas-mass-to-light ratio beyond r5oo. For the void and fila- 
ment directions, the IMLR profiles should increase with radius 
beyond r5oo. 

At r2oo, the cumulative IMLR was 0.007-0.01 Mq/Lk,q- 
A previously reported theoretical model predicts that the 
oxygen-mass-to-light ratio (OMLR) of a cluster is a sensi- 
tive function of the slope of the initial mass function (IMF) 
(Renzini 2005). Here, the oxygen mass is a sum of that 
trapped in stars and that in the ICM. By adopting a Salpeter 
IMF with a slope of 2.35, the expected value of the OMLR 
is ^ 0.1 Mq/Lb,q- In contrast, a top-heavy IMF with a 
slope of 1.35 overproduces metals more than that with a fac- 
tor of 20. For the Hydra A cluster, by adopting the ra- 
tio of stellar-mass-to-light ratio 4-6 in the B-band and K- 
band (Nagino & Matsushita 2009; Matsushita et al. 201 1), the 
IMLR in the B-band becomes 0.03-0.06 AIq/Lb.q and the 
cumulative OMLR within the virial radius becomes -^0. 1-0.3 



{Ao/Apc) Mq/Lb,q- Here, Aq and Apc are abundances of O 
and Fe in solar units, respectively. Beyond rsoo of the Hydra 
A cluster, the integrated stellar mass is only several percent of 
the integrated gas mass. Thus, even by considering the differ- 
ence of stellar metallicity and the abundance of the ICM, most 
of the metals in the cluster are in the ICM. Therefore, for cases 
with solar O/Fe ratio, the OMLR is similar to the expectations 
reported by Renzini (2005), assuming a Salpeter IMF in the 
Hydra A cluster In contrast, by adopting the O/Fe ratio from 
a nucleosynthesis model of SN II with metallicity = 0.02, as 
reported by Nomoto et al. (2006), the OMLR within the virial 
radius becomes 0.3-1.0 il/0/iB0 in the Hydra A cluster This 
value is larger than the expectation from the Salpeter IMF, and 
is consistent with a flatter IMF slope. 

In Figure 15, the cumulative IMLR profile of the Hydra A 
cluster was compared with those of other groups and clusters, 
including Coma (8 keV; Matsushita et al. 2011), Centaurus 
(4 keV; Sakuma et al. 2011), and AWM 7 clusters (3.6 keV; 
Sato et al. 2008), as well as NGC 1550 (1.2 keV; Sato et al. 
2010) and NGC 5044 groups (1.0 keV; Komiyama et al. 2009). 
Within 0.5r2oo, the cumulative IMLR profile of the Hydra A 
cluster increased with radius by a factor of two higher than for 
other systems. No systematic dependence on the ICM tem- 
perature was evident for clusters with temperatures higher than 
2 keV, although beyond 0.1r2oo, a group of galaxies shows a 
significantly smaller IMLR 

The stellar and gas-mass fractions within r^oo depend on the 
total system mass (Lin, Mohr, Stanford 2003; Lin & Mohr 
2004; Vikhlinin et al. 2006; Sun et al. 2009; Giodini et al. 
2009). These studies determined that within rsoo the stellar- 
to-total-mass ratios of the groups are much larger than those in 
the clusters, whereas the gas-mass fraction increases with the 
system mass. The observed higher stellar mass fraction and the 
lower gas mass fraction within r^oo in poor systems are occa- 
sionally interpreted as proof that the star formation efficiency 
depends on the system mass. However, as shown in Figure 
11, the gas is more extended than the stars in a cluster within 
rsoQ. The gas density profiles in the central regions of groups 
and poor clusters were observed to be shallower than those in 
the self-similar model, and the relative entropy level was corre- 
spondingly higher than that in rich clusters (Ponman, Cannon, 
Navarro 1999; Ponman, Sanderson, Finoguenov 2003; Sun 
et al. 2009). Then, the difference in the ratio of gas-mass- to- 
stellar-mass may reflect differences in distributions of gas and 
stars, which in turn reflects the history of energy injection from 
galaxies to the ICM. To study the fractions of stars and gas in 
clusters of galaxies, and its dependence on the system mass, 
we need measurements of gas and stellar mass beyond rsoo of 
other clusters. 

6. Summary and Conclusion 

We presented the results of the Suzaku observation of the 
Hydra A galaxy cluster, which is the first medium-sized clus- 
ter (temperature ~ 3 keV) observed with Suzaku out to the 
virial radius. Two observations were conducted, north-west 
and north-east offsets, which continue in a filament direction 
and a void direction of the large-scale structure of the Universe, 
respectively. We investigated possible azimuthal variations 
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in temperature, electron density, entropy and total mass pro- 
files. Our analysis revealed that distributions of X-ray emis- 
sion and galaxies elongate in the filament direction. The en- 
tropy profiles become flatter beyond r5oo in contrast to the 
r^'^ relationship expected from the shock heating model of 
the ICM. The entropy profiles are universal in clusters ob- 
served with Suzaku when scaled with the averaged tempera- 
ture of each cluster The hydrostatic masses in the two direc- 
tions agree and the NFW universal matter profile represents 
the hydrostatic mass distribution obtained up to 1.8 r5oo, with 
Cvi, = 9-10 and Mvir = (2.2 ± 0.4) x 10^'^Mq. The gas frac- 
tion, Mgas(< ?')/AfH.E.(< r) significantly exceeds the cosmic 
mean baryon fraction of WMAP7 beyond ^ t^qq. The flat- 
tening of the entropy profile and higher gas fraction contra- 
dict expectations based on numerical simulations. An under- 
estimate of gas temperature due to the discrepancy between 
ion and electron temperatures gives higher entropy and smaller 
gas fraction at cluster outskirts. However, a flat or increasing 
temperature profile is required. If bulk motions caused by in- 
falling matter from filaments of clusters are higher than that of 
numerical simulations, and if the ICM deviates from the hy- 
drostatic equilibrium, the entropy becomes smaller and the gas 
fraction is overestimated. In addition, we obtained IMLR up 
to the virial radius for the first time, and compared our results 
with other clusters and groups. The IMLR profile is consistent 
within 0.5 r2oo with other clusters, and become flatter from 0.5 
r2oo tor2oo- 
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Abstract 



o 
u 

Qh, We report Suzaku observations of the northern half of the Hydra A cluster out to ~ 1 .4 Mpc, reaching the virial 

radius. This is the first Suzaku observations of a medium-size {kT ~3 keV) cluster out to the virial radius. Two 
^ ' observations were conducted, north-west and north-east offsets, which continue in a filament direction and a void 

jXi ' direction of the large-scale structure of the Universe, respectively. The X-ray emission and distribution of galaxies 

C^ \ elongate in the filament direction. The temperature profiles in the two directions are mostly consistent with each 

other within the error bars and drop to 1 .5 keV at 1 .5 rsoo . As observed by Suzaku in hot clusters, the entropy profile 
tH- ' becomes flatter beyond rsoo, in disagreement with the r'^' relationship that is expected from accretion shock heating 

J> . models. When scaled with the average intracluster medium (ICM) temperature, the entropy profiles of clusters 

^D ' observed with Suzaku are universal and do not depend on system mass. The hydrostatic mass values in the void 

^^ \ and filament directions are in good agreement, and the Navarro, Frenk, and White universal mass profile represents 

the hydrostatic mass distribution up to ^ 2 7'5oo. Beyond rsoo, the ratio of gas mass to hydrostatic mass exceeds 
the result of the Wilkinson microwave anisotropy probe, and at rioo, these ratios in the filament and void directions 
C^ . reach 0.4 and 0.3, respectively. We discuss possible deviations from hydrostatic equilibrium at cluster outskirts. We 

derived radial profiles of the gas-mass-to-light ratio and iron-mass-to-light ratio out to the virial radius. Within rsoo, 
the iron-mass-to-light ratio of the Hydra A cluster was compared with those in other clusters observed with Suzaku. 
Key words: galaxies: clusters: individual (Hydra A cluster) — X-rays: galaxies: clusters — intergalactic 
medium 



.^.' 1. Introduction Imaging Spectrometer (XIS; Koyama et al. 2007), Suzaku 

(Mitsuda et al. 2007) was able to unveil for the first time 

Clusters of galaxies are the largest self-gravitating systems the intracluster medium (ICM) beyond rsoo, a radius within 

in the Universe, and offer unique information on the process of which the mean cluster-mass density is 500 times the cosmic 

structure formation governed by cold dark matter (CDM). In critical density. The accurate calibration of the XIS also al- 

addition, these clusters are considered as a laboratory for study- lows precise measurements of the ICM temperature (Sato et 

ing thermal and chemical evolutions of the Universe in which al. 201 1). Suzaku derived the temperature and entropy profiles 

baryons play the most important role. X-ray observations pro- of the ICM of several massive clusters up to the virial radius 

vide valuable information about the structure formation, gas (George et al. 2009; Reiprich et al. 2009; Bautz et al. 2009; 

heating and cooling, and metal enrichment of galaxy clusters. Kawaharada et al. 2010; Hoshino et al. 2010; Simionescu et 

Because the dynamical time-scale of clusters is comparable al. 2011; Akamatsu et al. 2011). From the center to r2oo, a 

to the Hubble time, cluster outskirts should maintain original systematic drop in temperature was found by a factor of ^ 

records of cluster evolution via accretion of gas and substruc- 3, and the observed entropy profiles become flatter beyond 

tures from filaments of the surrounding large-scale structure of rsoo. These profiles are lower than the entropy profile pre- 

the Universe. dieted by the numerical simulations of gravitational collapse 

Thanks to the low and stable background of the X-ray (Tozzi & Norman 2001; Voit, Kay, Bryan 2005), which is pro- 
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portional to r^'^. Kawaharada et al. (2010) discovered that be- 
yond rsoo of the Abell 1689 cluster, the total mass obtained 
from weak-lensing observations with Subaru is larger than that 
calculated assuming hydrostatic equilibrium. Therefore, one 
explanation for the low entropy profiles at cluster outskirts is 
that infalling matter retained some of its kinetic energy in bulk 
motion (Bautz et al. 2009; George et al. 2009; Kawaharada et 
al. 2010). Based on Suzaku observations of the Perseus cluster, 
Simionescu et al. (201 1) proposed a gas-clumping effect as an 
additional interpretation. Hoshino et al. (2010) and Akamatsu 
et al. (2011) discussed possible deviations of electron temper- 
ature from ion temperature to explain the observed lower tem- 
perature and entropy profiles. To clarify these effects in the out- 
skirts, dependence on the system mass or on the ICM temper- 
ature should be examined. With XMM-Newton observations. 
Urban et al. (2011) found a similar flattening of the entropy 
profile in one direction in the Virgo cluster {kT = 2.3 keV). 
However, Humphrey et al. (201 1) detected no evidence of flat 
profiles at large scales (> rsoo) in the relaxed fossil group or in 
the poor cluster, RXJ 1 159H-5531. Therefore, more samples of 
medium-sized clusters are required. 

In addition to measuring the temperature and entropy pro- 
files, Suzaku measured the abundance of Fe in the ICM beyond 
0.5ri8o (Fujita et al. 2008; Tawa et al. 2008; Simionescu et al. 
2011). Metal abundances in the ICM also provide important 
information on the chemical history and evolution of clusters. 
The ASCA satellite first measured the distribution of Fe in the 
ICM (Fukazawa et al. 2000; Finoguenov et al. 2001). Recently, 
spatial distributions of Fe from 0.3 to 0.4 riso have been stud- 
ied with the XMM-Newton and Chandra satellites (Vikhlinin 
et al. 2005; Maughan et al. 2008; Leccardi & Molendi 2008; 
Matsushita 201 1). Since metals are synthesized by supernovae 
(SNe) in galaxies, the ratios of metal mass in the ICM to the to- 
tal light from galaxies in clusters or groups, (i.e., metal-mass- 
to-light ratios) are key parameters in investigating the chemi- 
cal evolution of the ICM. Suzaku measured the iron-mass-to- 
light ratios (IMLR) of several clusters and galaxy groups out 
to 0.2 - 0.5 ri8o (Matsushita 2007; Komiyama et al. 2009; 
Sato et al. 2007; Sato et al. 2008; Sato et al. 2009a; Sato et al. 
2009b; Sato et al. 2010; Sakuma et al. 201 1), and with XMM- 
Newton, the IMLR of the Coma cluster out to 0.5ri8o was de- 
rived (Matsushita et al. 201 1). The IMLR profiles increase with 
radius, indicating that Fe in the ICM extends farther than stars. 

The Hydra A cluster (z = 0.0539) with an ICM tempera- 
ture of ^3 keV is one of the prototype cool-core clusters in 
which McNamara et al (2000) discovered a displacement of 
X-ray gas in the central region through the radio lobes from 
the central active galactic nucleus (AGN). This cluster is also 
known as Abell 780 and remains a major examples of AGN 
interaction that is studied through radio and X-rays (Taylor 
et al. 1990; David et al. 2001; Nulsen et al. 2002; Nulsen et 
al. 2005; Lane et al. 2004; Wise et al. 2007; Simionescu et 
al. 2009a; Simionescu et al. 2009b; Kirkpatrick et al. 2009). 
Through Chandra and XMM-Newton observations, a sharp 
X-ray surface-brightness edge was detected at radii between 
4.3' and 6' (200-300 kpc) and was interpreted as a shock 
wave caused by an AGN outburst. The abundances of the 
Hydra A cluster was measured within ~0.3 riso with XMM 
(Simionescu et al. 2009b; Matsushita 201 1). 
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Fig. 1. XISO image of the Hydra A cluster within 0.5^.0 keV energy 
range. The Ught blue cross indicates the X-ray peak of the cluster. 
Differences in the vignetting effect are uncorrected. Regions of spectral 
accumulation are shown as magenta rings. The ring radii are 2', 4', 6', 
8', 10', 12', 15' and 23'. Black circles represent excluded regions 
around point sources. 

This paper reports the results of two Suzaku observations 
in the northern half of the Hydra A cluster out to 25'(~1.5 
Mpc), which corresponds to the virial radius. The observa- 
tions were conducted using the XIS. This study reports the first 
Suzaku observations of a medium-sized cluster with an aver- 
age temperature of ^3 keV up to the virial radius. The two 
observed fields continue in a filament direction that continues 
to the Abell 754 cluster and a void direction of the large-scale 
structure of the Universe (see Figure 9). 

In this study, we use the Hubble constant, Hq = 
70 km s"^ Mpc^i, f7,„ = 0.27 and Qa = 0.73. The luminos- 
ity distance Z3l and angular size distance Da to the Hydra A 
cluster are I?l — 240 Mpc and Da = 217 Mpc, respectively, 
and 1' corresponds to 66.6 kpc. We use the epoch J2000.0 for 
the right ascension and declination of the equatorial coordi- 
nate system. In addition, we use the solar abundance table by 
Lodders (2003), in which the solar Fe abundance relative to H 
is 2.95 X 10^^. Considering a difference in solar He abundance, 
the Fe abundance yielded by Lodders (2003) is 1 .5 times higher 
than that using the photo-spheric value by Anders & Grevesse 
(1989). Errors were quoted at a 90% confidence level for a 
single parameter 

2. Observation and Data Reduction 

2.1. Suzaku observations 

Suzaku conducted two observations of the Hydra A cluster 
in November 2010, which was during the Suzaku Phase- V pe- 
riod. The details of the observations are summarized in Table 1 . 
The first observation, HYDRA A-1, was 12' north-west offset 
from the X-ray peak of the Hydra A cluster with coordinates 
(RA, Dec) = (139.5236, -12.0955) in degrees. The second ob- 
servation, HYDRA A-2, was 12' north-east offset from the X- 
ray peak of the Hydra A cluster The north-west and north-east 
fields continue into the filament and void structures, respec- 
tively (See figure 9). Hereafter, we refer to the north-west and 
north-east offsets as filament and void, respectively. 
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Suzaku Observations of the Hydra A Cluster 

Table 1. Observation log of the Hydra A cluster 



Field name 



Target name 



Sequence 
number 



Observation 
date 



RA* 

(deg) 



Dec* 

(deg) 



Exposure 

(ks) 



North-west (filament) 
North-east (void) 



HYDRA A- 1 
HYDRA A-2 



805007010 
805008010 



2010-11-08 
2010-11-09 



139.3723 
139.6749 



-11.9472 
-11.9472 



35.6 
34.9 



Coordinates are referred to J2000.0. 



The XIS, which was operated in its normal mode during the 
observations, consists of three sets of X-ray CCD (XISO, XIS 1, 
and XIS3). XISl is a back-illuminated (BI) sensor, while both 
XISO and XIS3 are front-illuminated (FI) sensors. Data re- 
duction was done with HEAsoft version 6.11. The XIS event 
lists created by rev 2.5 pipeline processing were filtered using 
the following additional criteria: a geomagnetic cutoff rigid- 
ity (COR2) >6 GV, and an elevation angle > 10° from the 
earth limb. The 5x5 and 3x3 editing modes data formats were 
added. The exposure times after data selection are shown in 
Table 1. 

To obtain the temperature and electron density profiles, we 
accumulated spectra within annular regions centered on the X- 
ray peak as shown in Figure 1. Regions around calibration 
sources and the innermost region within 2' were excluded from 
the spectral analysis. The reason for this exclusion is that we 
could not obtain sufficient photon statistics to analyze the spec- 
tra (which have multiphase features at the cluster center) be- 
cause the given area was detected on the edge of a CCD chip in 
our observation. In addition, Suzaku's PSF is also insufficient 
to resolve such a complex structure. However, with their pow- 
erful imaging capability, Chandra and XMM observations have 
unveiled such complex features in the central region. Using the 
ewavelet tool in SAS\ we searched for point-like sources in 
the Suzaku images in the energy range of 2.0-10.0 keV. In this 
energy range, the flux levels were approximately > 2 x 10^^^ 
erg cm^^ s^^. These point sources were excluded as circular 
regions with radii of 1' (figure 1). The non X-ray background 
(NXB) was subtracted from each spectrum using a database of 
night Earth observations with the same detector area and COR 
distribution (Tawa et al. 2008). 

We included the degradation of energy resolution due to ra- 
diation damage in the redistribution matrix file (RMF) gener- 
ated by the xisrmfgen Ftools task. In addition, we created 
an ancillary response file (ARE) using the xissimarfgen 
Ftools task (Ishisaki et al. 2007). A decrease in the low-energy 
transmission of the XIS optical blocking filter (ORE) was also 
included in the ARE. For filament and void fields, we generated 
ARE files assuming uniformly extended emission from an en- 
riched region with a 20' radius. We used the XSPEC_vl 2.7.0 
package and ATOMDB_v2.0.1 for spectral analysis. Each 
spectrum was binned and each spectral bin contained a min- 
imum of 50 counts. To avoid systematic uncertainties in the 
background, we ignored energy ranges above 7 keV and below 
0.7 keV. In addition, we excluded the narrow energy band be- 
tween 1.82 and 1.84 keV in the fits because of the incomplete 
response around the Si edge. 
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2.2. XMM observations 

We used the same data and analysis as that reported by 
Matsushita (201 1). The XMM-Newton archival data of Hydra 
A (observation identifier 0109980301) had exposure times of 
14.0 ks, 19.1 ks and 22.2 ks forMOSl, MOS2 and PN, respec- 
tively, after background flares were screened out. We selected 
events with patterns smaller than 5 and 1 3 for the PN and MOS, 
respectively. Spectra were accumulated in concentric annular 
regions of 0'-0.5', 0.5'-l', l'-1.5', 1.5'-2', 2'-3', 3'-4', 4'-5', 
5'-6', 6'-7', 7'-8', and 8'-9', centered on the X-ray peak of the 
Hydra A cluster Here, the X-ray peaks were derived using the 
ewavelet tool of SAS-v8.0.0, and luminous point sources were 
excluded. The spectra from MOSl and MOS2 were added. 
The background spectrum for each annular region was calcu- 
lated by integrating blank-sky data in the same detector region. 
From deep-sky observations with XMM-Newton, we selected 
the data having background most similar to that of the Hydra 
A cluster and the faintest Galactic emission, after screening 
out background flare events from the data and the background, 
following Katayama et al. (2004). Next, we scaled the back- 
ground using a count rate between 10 and 12 keV. 

The response matrix file and the ARE corresponding to each 
spectrum were calculated using SAS-v8.0.0. Further details 
appear in sections 2 and 3 of Matsushita (201 1). 

3. Data Analysis 

3.1. Estimation of background spectra 

To derive the radial profiles of the temperature, electron 
density, and Ee abundance of the ICM, we fitted the NXB- 
subtracted spectra with a thermal plasma model (APEC: Smith 
etal. 2001). 

We first fitted the spectra in the outermost ring (15'-23' re- 
gion) to determine the local X-ray background. As reported 
by Yoshino et al. (2009), the background emission of Suzaku 
XIS can be fitted with a three-component model: two thermal 
plasma models (APEC; Smith et al. 2001) for the local hot 
bubble (LHB) and the solar wind charge exchange (SWCX), 
the Milky Way Halo (MWH), and a power-law model for 
the extragalactic cosmic X-ray background (CXB). MWH and 
CXB components were convolved with an photoelectric ab- 
sorption in the Galaxy, A^h ■ Therefore, we adopted the model 
formula, apecLHB + wabs x (apecMWH + power-lawcxB) as 
background. We assumed a zero redshift and a solar abundance 
for LHB and MWH components. The temperature of the LHB 
and MWH was fixed at 0.1 keV and 0.3 keV, respectively. The 
column density of the Galactic neutral hydrogen was fixed at 
4.7x1020 cm"2 (Kalberla et al. 2005), and the photon index 
of the CXB component was fixed at 1.41 (Kushino et al. 2002). 
Using the formula apecLHB + wabs x {apecicM + apecuwn + 
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Table 2. Resultant parameters of the background components. 



Cosmic X-ray Background Local Hot Bubble MiUcy-Way halo 



normalization* 



normalization^ 



normalization^^ 

xlO-4 



.92: 



-0.85 



.93: 



-5.75 



1.83: 



-1.72 



Reduced-x'^ 

X^/d.o.f 



1.04(144/139) 



* Measured at IkeV with unit of photons cm^'^ s^'^ keV^'^ sr^'^. 

t Normalization of the APEC, component divided by the solid angle, f!'^, assumed in the uniform-sky ARF calculation (20' radius), Norm 
f 71^11 jjdV/{4:TT{l + z)'^D^)/Q" X 10"^'* cm"'' 4007r arcmin^^, where D^j is the angular distance to the source. 



power-lawcxB), we fitted the spectra with a sum of the back- 
ground components and a single-temperature APEC model for 
the ICM. Temperature and normalization of the ICM compo- 
nent were allowed to vary, and redshift was fixed at 0.0539. 
Normalizations of the background components were also left 
free. The spectra of the filament and void fields were fitted 
simultaneously, and each parameter of background and abun- 
dance of the ICM was assumed to be the same. 

Results of the spectral fit are shown in Figure 2. The result- 
ing parameters for the background and ICM are shown in Table 
2 and 3, respectively. The normalizations obtained for the CXB 
and the Galactic components are consistent with those derived 
by Kushino et al. (2002) and Yoshino et al. (2009), respectively. 

Several blank fields observed with Suzaku contained emis- 
sion with kT = 0.6 - 0.8 keV (Yoshino et al. 2009). Therefore, 
we added an APEC component with kT = 0.6 keV and refit- 
ted the outermost spectra. The normalization of the 0.6 keV 
component had a much smaller value than the other Galactic 
components, and the temperatures and normalizations of the 
ICM did not change. 
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Fig. 2. NXB-subtracted spectra of XIS 1 (crosses) for 
the outermost ring (15'-23' region), fitted by the 
apecLHB + wabs X (apecjcM + "-pecMwu + power-lawcys.) 
model (stepped solid lines). Black and red colors con'cspond to the 
directions in the filament and void, respectively. The bottom panel 
shows residuals of the fit. The contributions of the ICM are plotted as 
black and red solid lines, those of the Galactic emission as orange and 
green fines, and those of the CXB as blue lines. 
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Fig. 4. Radial temperature profiles of the filament (black) and 
void (red) directions. Solid lines show the best fit function of 
A X exp(— fe X r). Chandra results (David et al. 2001) are plotted in 
blue. Vertical dotted fines show rsooi^iso. and rioo, derived from 
the fitting of the hydrostatic mass with the NFW universal mass profile 
(Navarro et al. 1997), as described in subsection 4.4. 

3.2. Spectral fits of annular regions 

We fitted the XIS spectra in each region of the Hydra A clus- 
ter using a single-temperature APEC model for the ICM with 
the Galactic absorption, A^h- The spectra of each annular re- 
gion of the two fields were fitted simultaneously by minimiz- 
ing the total x^ value. Here, the common background model 
was included for all regions, where the surface brightness of 
the background components were restricted within the statis- 
tical errors derived for the outermost regions. The tempera- 
ture and normalization of the ICM component were free pa- 
rameters. The abundances in each annular region in the two 
fields were assumed to have the same value. A^h was fixed to 
the Galactic value of 4.7 x 10'^° cm^^ in the direction of the 
Hydra A cluster 

The resultant parameters are summarized in Table 3, and the 
best-fit spectra are shown in Figure 3. The spectra are well 
represented with the single-temperature model for the ICM and 
the background. 

4. Results 

4.1. Temperature profile 

Figure 4 shows the radial temperature profiles of the two di- 
rections derived from the spectral fits. At a given radius, the 
temperatures for each direction are consistent within statisti- 
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Fig. 3. NXB-subtracted XISl spectra fitted with the single temperature (APEC) model for the ICM. Black and red colors correspond to directions toward 
the filament and void, respectively. Lower panels in each diagram show residual of fit. The contributions of the ICM are plotted as red and black solid 
Unes, while those of the galactic emission and the CXB are plotted as orange, green and blue lines, respectively. 
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Table 3. Parameters resulting from fitting ICM with the single-temperature APEC model in the energy ranges of 0.7-7.0 keV. 



[Vol. , 



Region 



filament 



void 



(arcmin) 



(kpc) 



2' -4' 133 - 266 3.81+^;^ 



kT (keV) abundance (solar) normalization* kT (keV) normalization* 



-0.02 

0^^0+0.01 



-0.14 

3.09+°-2^ 



-0.02 

0.0901H1 
0.035+2:^3 

nni5+°''03 

U.U10_o 0Q2 

fO.0016 

-0.0012 

S+0.0012 

-0.0012 

002^+''°°" 



X /d.o.f. 



4' -6' 
6' -8' 
8' - 10' 
10' - 12' 
12' - 15' 
15' - 23' 



266 - 400 
400 - 533 
533 - 666 
666 - 799 
799 - 997 
997 - 1529 



-0.34 
q ni+0.32 
'J"^^-0.23 
q 11+0.28 
'-'■^^-0.23 

Z.OU_(j 27 

2 '^.A+O-a^ 

^■'J^-0.36 

qf? + 0.33 

1 44+0.47 
^■^^-0.33 



0.33 ±0.10 
0.43 ±0.43 
0.28 ±0.15 
0.18 ±0.13 
0.38 ±0.38 
0.30 ±0.19 
0.10 ±0.08 



0.45-* 



3.18^ 



-0.01 

U.UOO_Q (jog 

U.UZO_Q QQ3 

r,Qio+0.003 

010+°°''^ 

U.UlU_o oQj^ 

0036+''°°^2 



-0.24 
7+0.38 

-0.31 
:;+0.86 

-0.81 

4(5+1.60 

"1+0.33 
-0.31 

1 Qq+0.65 
-'-■•J'^-0.42 



3.07! 
3.161 



1.79! 



0.36-* 



0.0058: 
0.0048^ 



1.22(263/216) 
1.09(146/134) 
1.18(177/150) 
1.18(129/109) 
1.07(99/93) 
1.13(170/150) 
1.19(263/221) 



* Nomiahzation of the APEC, component divided by the solid angle, Q", assumed in the unifomi-sky ARF calculation (20' radius), Norm ■ 
r?ie?i£fdy/(47r(l + z)^D'^)/fl" X 10^^* cm"''' 4007r arcmin"^, where D^ is the angular distance to the source. 



cal errors, except around ^200 kpc from the center The re- 
gion around 200 kpc corresponds to the shock front detected 
by Chandra and XMM in the filament direction (Simionescu et 
al. 2009b; Gitti et al. 2011). From 300 kpc to 1200 kpc, the 
temperatures in each direction decrease with radius in a similar 
manner down to ^ 1 .4 ke V. 

From ROSAT observations, Ikebe et al. (1997) found that the 
spectrum at the central region within 2' exhibits the multiphase 
feature. Although the central 2' region was excluded in our 
spectral analysis, we examined the spectral fits with the two- 
temperature model in regions with r > 2'. The reduced-^^ 
for the two-temperature model was not significantly improved 
in comparison with the single-temperature model. 

To determine the hydrostatic mass of Hydra A, we fitted the 
temperature profile of each direction with the exponential for- 
mula, A X exp(— 6 X r), where A and b are free parameters with 
units of keV and kpc"^, respectively, and r is the distance from 
the X-ray peak of the cluster center The resulting fit parame- 
ters and best-fit functions are shown in Figure 4. In the filament 
direction, the normaUzation factor A of the exponential model 
was 18% higher than that in the void direction, which reflects a 
temperature difference of ^ 200 kpc. The parameter b for both 
directions is consistent within statistical error 

4.2. Deprojected electron density profile 

We calculated three-dimensional electron density profiles 
from the normalization of the ICM component derived with 
Suzaku, considering the geometrical volume that contributes 
to each two-dimensional region. Here, we assumed spherical 
symmetry within each field. The electron density profile was 
also obtained using XMM data. To avoid uncertainties in the 
background, we fitted the annular spectra of XMM MOS with 
the APEC model in an energy range of 1.6-5.0 keV. Here, the 
temperature and Fe abundance were left free within the error 
bars derived by Matsushita (201 1). We fitted the derived emis- 
sivity profiles with a double-/3-model, and calculated the three- 
dimensional radial profile of the electron-density. 

Figure 5 shows the electron density profile derived from 
XMM and Suzaku (this work), Chandra (David et al. 2001) and 
the ROSAT surface-brightness profile scaled with the XMM 
density profile at 200 kpc (Ikebe et al. 1997). The deprojected 
electron density profile derived from XMM and Suzaku obser- 
vations agree well with those from ROSAT within 500 kpc and 
from Chandra within 200 kpc. 

At a given distance from the X-ray peak, electron density 



LU 




100 200 
distance (kpc) 

Fig. 5. Deprojected electron density profile of the filament (black) and 
void (red) directions. Black and red solid lines show the best fit with 
the double-/? model. Green solid line shows the electron density profile 
derived with XMM data, the light blue solid line shows that derived 
from the ROSAT surface-brightness profile (Ikebe et al. 1997) scaled 
with the XMM profile, and the blue filled circles show that derived 
from Chandra results (David et al. 2001). 

beyond 400 kpc in the filament direction was systematically 
higher than that in the void direction, reflecting the differ- 
ence in normalizations (Table 3 and Figure 3). At 900 kpc, 
for example, the deprojected electron densities in the fila- 
ment and void directions were (2.54 ± 0.18) x 10~* cm~^ and 
(1.76 ± 0.21) X lO""* cm"^, respectively. The electron density 
dropped to ^lO""' cm~^ around 1200 kpc. 

We fitted the density profiles from XMM and Suzaku results 
with the double-/? model from 40 to 1200 kpc. The results 
are shown in Figure 5, and the resulting parameters are given 
in Table 4. From 400 kpc to 1530 kpc, the radial profiles of 
the electron density are also well represented by a power-law 
profile with slopes of r = 1.60 ±0.05 and 1.72±0.07in the fil- 
ament and void directions, respectively. The slope of the elec- 
tron density in the void direction was slightly steeper than that 
in the filament direction. Fitting the electron density profiles 
beyond 100 kpc with the single /3-model resulted in /3 = 0.61 
± 0.02 and (i = 0.1 A ± 0.03 for the filament and void direc- 
tions, respectively. These values are close to the value of 0.68 
derived from the ROSAT surface-brightness profile (Ikebe et 
al. 1997). 



No. ] Suzaku Observations of the Hydra A Cluster 

Table 4. Parameters resulting from fitting tlie deprojected density profile derived from XMM and Suzaku with the double-/3 model*. 



Field 



narrower component 



wider component 



710,1 (10 cm ■=') rc,i (kpc) 



/3i 



no,2 (10 -^cm ■^) rc.2 (kpc) 



/32 



void 
filament 



6.73 ±0.17 
6.61 ±0.19 



21.3±4.7 
22.3±5.4 



0.59 ±0.01 
0.62 ±0.01 



2.16±0.14 
2.73±0.15 



278 ± 264 
193 ±114 



0.95 ±0.07 
0.60 ±0.03 



* double-/3 model is given by the form n{r) = no,i(l + (r/rc^i)'^) '^'^i/'^ + no, 2(1 + (''"/f'c,2)'^) 
The errors of no,i and /3i are calculated by fixing the rc,i and wider component, and vice versa. 



-3^/2 



Table 5. Resultant fit pai'ameters of the deprojected density profile derived from XMM and Suzaku with the single-/? model* beyond 100 kpc. 



Field no (lO"'' cm"^) 



(kpc) 



(3 



void 7.38 ±0.42 160 ±13 0.74 ±0.03 

filament 7.98 ±0.50 131 ± 10 0.61 ±0.02 

' single-/? model is given by the form n{r) = no(l + {r j r c)^')^'^^ ^ '^ . 
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Fig. 6. Entropy profiles in the filament (black) and void (red) direc- 
tions. Crosses and diamonds show entropy calculated from the result- 
ing temperature and electron density, respectively, as shown in Figures 
4 and 5. Chandra results (David et al. 2001) are represented by blue 
crosses. The green line shows the entropy profile calculated by combin- 
ing the electron density from XMM and the temperature from Suzaku. 
The light blue line extrapolates the Suzaku and XMM data (green line) 
with a power-law formula fit for the data beyond 200 kpc, which has a 
fixed index of 1 . 1 . 

4. 3. Entropy profiles 

Entropy, which is a useful tool for investigating the thermo- 
dynamic history of hot gas in groups and clusters, is defined 
as 

kT 



S^ 



2/3' 



(1) 



where T and rie are the temperature and deprojected electron 
density, respectively. In Figure 6, we show the entropy pro- 
files of the Hydra A cluster calculated from the derived tem- 
perature and electron density profiles shown in Figures 4 and 
5. Entropy in the two directions increased from the center to 
800 kpc, or rsoo, which was derived from the Navarro, Frenk 
ans White(NFW) model fit to the hydrostatic mass described in 
section 4.4. Between 500 kpc and 800 kpc, the entropy in the 
void direction was systematically higher than in the filament 



direction. At the two outermost radial bins (beyond rsoo), the 
entropies of the two directions agree well. In addition, we cal- 
culated the entropy profile by combining the parameters result- 
ing from the fit of the electron density profile from XMM and 
the temperature profile from Suzaku as indicated by the green 
line in Figure 6. The entropy profiles using Suzaku, XMM, and 
Chandra data are consistent with each other (figure 6). 

The derived entropy profiles were compared with a power- 
law model with a fixed index of 1.1, which was expected from 
simulations of the accretion shock heating model (Tozzi & 
Norman 2001; Ponman, Sanderson, Finoguenov 2003). From 
200 kpc to 500 kpc, the derived entropy profiles agreed well 
with the r^^ relationship. In contrast, beyond 800 kpc, or 
^ rsoo, the entropy profiles in the two directions were system- 
atically lower than the predicted by the r^^ relationship. Inside 
200 kpc, the entropy exceeded the model prediction, which at- 
tribute to the presence of the cool core boundary (Peres et al. 
1998) and the shock front around 200 kpc. 

4.4. Hydrostatic mass 

We estimated the hydrostatic or gravitational mass assuming 
the hydrostatic equilibrium given by, 



MH.E.(<r) 



kT{r)r /dlnp{r) dlnkT{r) 



dlnr 



din? 



(2) 



where A/h.e.(< r) is the hydrostatic mass within the three- 
dimensional radius r, G is the gravitational constant, k is the 
Boltzmann constant, p is the gas density, kT is the tempera- 
ture, and yUJTip is the mean particle mass of the gas (the mean 
molecular weight is /i = 0.62). Here, we assumed spherically 
symmetric mass distribution. 

On the basis of the temperature and electron density pro- 
files, we calculated the hydrostatic mass of the Hydra A clus- 
ter beyond 100 kpc. The temperature gradient of the ith shell 
was calculated from a power-law fit of temperatures of shells 
(i — 1), i, and [i + 1). We employed the electron density pro- 
file from the single-/? model fit shown in table 5. The derived 
hydrostatic mass is plotted in Figure 7 as diamonds with the 
integrated gas mass profiles and stellar mass profiles derived in 
subsection 4.6. The derived mass was calculated for the whole 
azimuthal angle under the assumption of spherical symmetry 
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Fig. 7. left:Integrated hydrostatic and gas-mass profiles. All the black and red colors show the filament and void directions, respectively. The masses 
were calculated in the whole azimuthal angle under the assumption of spherical symmetry for each direction. The hydrostatic mass was calculated from 
the temperatures of shells (« — 1), i, and (i + 1) (diamonds) and the exponential relationship (dashed lines). All dashed lines show 90% upper and lower 
Umits. The green line shows the best-fit NFW model. The integrated stellar mass profiles in the filament and void directions using 2MASS data are 
represented by stars. Blue stars are the integrated stellar mass profile in the whole azimuthal angle coiTected for the contribution of the fainter galaxies 
using LF (solid). right:Cumulative profiles of gas fraction, ,fgas(<'') assuming hydrostatic equilibrium in both directions. The diamonds show /gas(< ^) 
of the integrated hydrostatic and gas mass of the corresponding directions. Dashed lines represent /gas (< r) applying the result of the NFW model fit. 
The value of the cosmic baryon fraction (light blue) from WMAP (Komatsu et al. 2011) is plotted. The stars represent the stellar fraction, /star(< t), 
using the NFW model fit. The meanings of the colors and fine styles are the same as in the left panel. 



using the parameters derived for each direction: thus, the mass 
in each direction was four times larger than that in the corre- 
sponding direction. Although the electron density in the fil- 
ament direction at a given radius was higher than that in the 
void direction, the difference in the slopes of the electron den- 
sity was small. As a result, the hydrostatic masses derived for 
the two directions agree well. In addition, we calculated the 
90% range of hydrostatic mass using the exponential relation- 
ship of the temperature profiles and fits of the electron-density 
profiles with the double-/? model, which are plotted as dashed 
lines in Figure 7. 

We fitted the derived hydrostatic mass profiles with the fol- 
lowing NFW formula, which shows the equilibrium density 
profile of the dark matter halo (Navarro et al. 1996; Navarro 
et al. 1997); 



5c = 



200 



(4) 



MNFw(<r)=47r(5cPcrs 



ln(l + a;) 



1 



,x = r/rs,(3) 



where AfNFw(< ^) is the mass within a radius r, p^ is the 
critical baryon density of the Universe, r^ is the scaled radius, 
and (5c is the characteristic density, that can be expressed in 
terms of the concentration parameter (c = tiqqJTs) as 



3 ln(l + c)-c/l + c' 

Since the hydrostatic masses derived for the two directions 
agree well, we fitted the two hydrostatic mass profiles with 
Suzaku and that with Chandra simultaneously, excluding the 
shock region from 100 kpc to 200 kpc. The parameters derived 
from the fit are shown in Table 7. The quantity rioo is equal 
to Ttjir in Hydra A cluster, since for our fiducial cosmologi- 
cal model, A„ir — 100 for halos with a redshift z ~ 0.0539 
(Nakamura & Suto 1997). We estimated r^ for various over- 
densities: A = 500, 200, 180 and 100, and ta derived from 
fits with the NFW model are 811 ± 10 kpc, 1189 ± 96 kpc, 
1243 ± 16 kpc, and 1577 ± 22 kpc, respectively. The rioo is 
close to the outermost ring radius. Possible deviations from hy- 
drostatic equilibrium have been previously discussed to explain 
the entropy flattening beyond r5oo (Bautz et al. 2009; George 
et al. 2009; Kawaharada et al. 2010): therefore, we fitted the 
hydrostatic mass out to tsoq. We were able to obtain almost 
the same results as shown in Table 7. Thus, the NFW model 
represents the hydrostatic mass out to rioo within error bars as 
shown in Figure 7. 
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Fig. 8. H.E. mass density profiles for ttie filament (black) and void 
(red) directions. Density was calculated from the integrated H.E. mass 
profile. The r~^ and r~^ profiles are plotted with solid blue and or- 
ange lines, respectively. 

Moreover, we calculated the total mass density profile from 
the total integrated hydrostatic mass profile in each direction. 
The density profile was also calculated from the NFW model 
obtained above. The calculated density profiles were compared 
with r^^ and r^^ profiles (Figure 8). The NFW density profile 
varied from r^^ to r^^ with increasing radius. Within r^Qo, 
the densities are well represented by the r^^ profile, which is 
shallower than the asymptotic matter density slope p ex r~^. 
Beyond rsoo, the r^^ profile is preferable than the r^^ profile. 

4.5. Gas fraction 

We derived the cumulative gas mass fraction as. 



/gas(<r) 



Mgas(<r) 

A/H.E.(<r)^ 



(5) 



where A/gas(< r) and A/h.e(< r) are the gas mass and hy- 
drostatic mass, respectively, within a sphere of radius r. As 
shown in Figure 7 and Table 6, the gas fraction increased in 
the outskirts and reached the cosmic mean baryon fraction de- 
rived from seven-year data of Wilkinson microwave anisotropy 
probe (WMAP7; Komatsu et al. 2011) at -- 0.5 tsoq. At 
'^looi /gas(< ?') in the filament direction exceeded the result 
of WMAP7, while the lower limit of /gas(< r) in the void 
direction reached the WMAP 7 fraction. Adopting the fit of 
the NFW model to the hydrostatic mass resulted in /gas(< r) 
in the void and the filament directions at rioo, exceeding the 
cosmic baryon fraction by a factor of 2 and 3, respectively. 
Furthermore, /gas(< r) in the void direction at rigg was close 
to that derived for the north-west direction of the Perseus clus- 
ter at r2oo (Simionescu et al. 2011): rigo was derived from the 
NFW model fit of the hydrostatic mass. The radial profile of 
/gas(< r) of the Abell 1246 cluster with kT = 6.0 keV exhib- 
ited similar behavior as that of the our result of the Hydra A 
cluster (Sato et al. 2012). 

4.6. K-band Luminosity of galaxies and stellar fraction 

The ratio of gas mass to stellar luminosity is a key pa- 
rameter for studying the star formation efficiency in clusters 



of galaxies. Since almost all metals in the ICM are synthe- 
sized in galaxies, the metal-mass-to-hght ratio provides a use- 
ful measure for studying the ICM chemical evolution. To es- 
timate the gas-mass and Fe-mass-to-light ratios, we calculated 
the A' -band luminosity of Hydra A cluster on the basis of the 
Two Micron All Sky Survey (2MASS) catalog, whereby all 
data is presented in a 2°x2° box: the coordinates of objects 
are shown in the left panel of Figure 9. The smoothed dis- 
tribution of galaxies around Hydra A obtained from 2MASS 
is shown in the right panel in Figure 9. The distribution of 
galaxies is elongated in the filament direction toward the Abell 
754 cluster, which exhibited a redshift close to that of Hydra 
A. In the void direction, the density of galaxies was signif- 
icantly lower. Toward the south, an additional filament-like 
distribution of galaxies appeared from the Hydra A cluster 
Hydra A itself was discovered to have an apparent magnitude 
of rriK = 11.07, or Lk = 5 x lO^^iif^© using a luminosity dis- 
tance Dl ~ 240 Mpc and a foreground Galactic extinction of 
Ak = 0.015 (Schlegel, Finkbeiner, Davis 1998). Among galax- 
ies in the NED data base^ with available redshifts, those with 
heliocentric velocities differ from that of the Hydra A cluster 
by >2000 km/s were excluded from the analysis. The /v -band 
surface brightness profile centered on Hydra A is shown in fig- 
ure 10. Beyond ^2000 kpc from the cluster center, the AT-band 
surface brightness flattened. Therefore, we adopted the region 
between 2000 kpc and 4000 kpc as a background. 

Figure 10 also shows the luminosity function (LF) of galax- 
ies in the Hydra A cluster within rioo- Here, we used as back- 
ground the LF of galaxies from 2000 kpc to 4000 kpc from the 
center of the Hydra A cluster The completeness limit of galax- 
ies of 2MASS is Kg = 13.5 in apparent magnitude ^, which 
corresponds to the i^-band luminosity, Lk=7 x lO^^L^.o for 
the distance of the Hydra A cluster The total A'-band lumi- 
nosity of galaxies above the 2MASS complete limit and that 
of all the detected galaxies including those below the limit are 
1.5 xIO^^Lk^o and 2.5 x 10^^ A^.o, respectively. To calcu- 
late the total light of the cluster, we evaluate the contribution 
of fainter galaxies below the 2MASS limit. The LF of galaxies 
in clusters can be described with the Schechter function as. 



4){L)dL = (f)^ 



t) -p(-vi.)<^ 



(6) 



where the parameters L», 0,^, and a are the characteristic lu- 
minosity, number densities, and faint-end power-law index, re- 
spectively. Upon adopting the A* of stacked cluster galaxies 
by Lin, Mohr, Stanford (2003) and integrating the Schechter 
function below the 2MASS limit, the total A'-band luminos- 
ity becomes 2.2 xlO^^A^.o, which is smaller than that of 
all the galaxies detected by 2MASS including those below 
the complete limit. If we allow A, to vary, the fit for the 
LF of the Hydra A cluster improves, as shown in Figure 10, 
and the total K-band luminosity becomes 2.8 xlO^^Ai^-©. 
Therefore, the simple sum of AT-band luminosities of the de- 
tected galaxies, including those below the 2MASS complete 
limit is close to that integrated luminosity of the Schechter 
function. Therefore, we adopted the total AT -band luminosity 



^ http://nedwww.ipac.caltech.edu/ 

^ See http://www.ipac.caltech.edu/2mass/releases/second/doc/explsup.html 
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Table 6. Cumulative values of the hydrostatic mass (Mh.e.), gas mass (Mgaa), and gas fraction (fgas) of the Hydra Hydra A cluster. 



[Vol. , 



Region 




filament 






void 




r 
(kpc) 


Mgas 

(10'-* Mq) 


Mh.e. 
(10^* Mq) 


Jgas 


Mgaa 
(10" Mq) 


Mh.e. 

(10" Mq) 


J gas 


133 - 266 


0.052 ±0.001 


0.39 ±0.04 


0.13 ±0.01 


0.052 ±0.001 


0.36 ±0.02 


0.15 ±0.01 


266 - 400 


0.11 ±0.01 


0.72 ±0.07 


0.15 ±0.02 


0.11 ±0.01 


0.68 ±0.06 


0.16 ±0.02 


400 - 533 
533 - 666 


0.19 ±0.02 
0.26 ±0.03 


1.06 ±0.09 
1.27 ±0.21 


0.17±0.02 
0.21 ±0.04 


0.17 ±0.02 
0.24 ±0.03 


1.03 ±0.13 
1.43 ±0.39 


0.17 ±0.03 
0.17 ±0.05 


666 - 799 
799 - 997 


0.35 ±0.05 
0.46 ±0.07 


1.35 ±0.35 
2.04 ±0.37 


0.26 ±0.08 
0.22 ±0.05 


0.30 ±0.04 
0.38 ±0.06 


3.03 ±1.41 
1.66 ±0.31 


0.10 ±0.05 
0.23 ±0.06 


997-1529 


0.71 ±0.12 


1.82 ±0.60 


0.39 ±0.14 


0.55 ±0.11 


1.67 ±0.78 


0.33 ±0.17 



Table 7. Resulting radii Tc and r2oo. and concentration parameter C200 derived from the NFW model. 



(kpc) r2oo (kpc) 



fitted up to rsoo 
fitted up to rioo 



126 ±8 
123 ±8 



1189 ±! 
1183 ±< 



9.44 ±0.48 
9.62 ±0.49 



of (2.5 ±0.3) xlQi^i^Q. 

We collected A' -band luminosity of galaxies in the north- 
west and north-east sectors for the filament and void direc- 
tions, respectively, as well as in the whole azimuthal angle and 
derived brightness profiles in the /f-band. For the whole az- 
imuthal angle, we collected galaxies above the 2MASS com- 
plete limit and evaluated the contribution of fainter galaxies 
by integrating the Schechter function. We obtained almost 
the same radial profile of A'-band luminosity as that from the 
simple sum of luminosities of all the galaxies detected with 
2MASS. Therefore, for the filament and void directions, we 
collected all the galaxies detected by 2MASS in each direction 
and derived radial profile of K-band luminosity. After subtract- 
ing the background, we deprojected the brightness profile of 
galaxies (under the assumption of spherical symmetry) and de- 
rived three-dimensional radial profiles of A'-band luminosity 
in the filament, void, and whole azimuthal angle. 

Figure 7 shows the integrated AT-band luminosity profiles. 
As opposed to the integrated luminosity of the whole azimuthal 
angle, figure 7 uses the four-times luminosity of the filament 
and void directions. The filament direction shows a system- 
atically higher A'-band luminosity than that in void direction. 
The integrated A' -band luminosity profiles of these two direc- 
tions become flatter from ~^800-1000 kpc, and the contribution 
of galaxies beyond 1000 kpc is small. The integrated luminos- 
ity profile of the whole azimuthal angle is between those in the 
filament and void directions except for the outermost region 
around 1200 kpc (Figure 7 and 9). This result is attributed to 
the luminosity of the whole azimuthal angle being affected by 
the filamentary structure around 1200 kpc in the south of the 
Hydra-A cluster 

We derived cumulative stellar mass fraction as. 



/star(< r) 



Mstar(< r) 



(7) 



MH.E.(<r) 

where, Mstar(< f) is the total stellar mass within r, derived 
from the total K-band luminosity within r, under the assump- 
tion that the stellar-mass-to-light ratio in the K-band is unity 
(Nagino & Matsushita 2009). At a given radius, the /star(< r) 



in the filament direction is higher by a factor of 2-3 than that 
in the void direction, reflecting the concentration of galaxies in 
the filament direction (Figure 7). Beyond r5oo, ./star(< r) is 
smaller than ^ 1%, and does not contribute to the baryon frac- 
tion. We also derived the stellar mass fraction using the whole 
azimuthal angle using the NFW mass for the filament and void 
directions. Although the stellar fraction continues to decrease 
with radius in the filament and the void directions, the stel- 
lar fraction in whole azimuthal directions shows a minimum at 
~ ^500, and increases with radius by a factor of '^ 1 .5 from rsoo 
tor 100. 

Using the integrated AT-band luminosity profile, we calcu- 
lated the gas-mass-to-light ratio in the filament and void direc- 
tions (Figure 11). The radial profile of the gas-mass-to-light 
ratio in the two directions increased with radius out to rioo. 
We also calculated the gas-mass-to-light ratio in the whole az- 
imuthal angle, using the stellar luminosity in the whole az- 
imuthal angle and the gas mass in the void direction which 
agrees that for the whole azimuthal angle derived by ROSAT 
(Ikebe et al. 1997). 

5. Discussions 

The Hydra A cluster is the first medium-sized cluster ob- 
served with Suzaku out to the virial radius. Suzaku observed 
the northern half of the cluster, and we derived radial profiles 
of temperature, electron density, entropy, hydrostatic mass, and 
K-band stellar luminosity in the filament and void directions. 
Temperature and entropy carry important information about 
the thermal history of the ICM. In subsection 5.1, we com- 
pare these profiles of the Hydra A cluster and other clusters 
observed with Suzaku and XMM to study their dependence 
on the system mass or on average ICM temperature. Next, 
we discuss possible explanations for entropy flattening and the 
higher gas-mass-to-hydrostatic-mass ratio at cluster outskirts. 
Discrepancies between electron and ion temperatures are dis- 
cussed in subsection 5.2, deviations from hydrostatic equilib- 
rium and the gas-clumping effect are discussed in subsections 
5.3 and 5.4. We also derive IMLR of the Hydra A cluster out 
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Fig. 9. Two-dimensional density distributions of galaxies around Hydra A, retrieved from 2MASS in the ftT-band. left: Stars represent galaxies with 
available redshift in the NASA/IPAC Extragalactic Database (NED). Green and blue stars indicate members and non-members, respectively, of the Hydra 
A cluster based on their redshift. The boxes represent the Suzaku fields of view (FOVs) of XIS pointings. The directions of the filament (north-west) and 
the void (north-east) are separated as dashed lines. Right: The box size is 5° x5°. Abell 754 cluster, with coordinates (RA, Dec) = (137.2087, -9.6366) 
in degrees, is shown in top right. A Gaussian smoothing scale is represented by the white hatched circle at bottom right. 
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Fig. 10. left:Surface-brightness profile of Hydra A observed with 2MASS photometric data in the A'-band. Blue points denote the detected galaxies of 
Hydra A and red stars denote the background-subtracted data. The dashed Une shows the subtracted background level and the dotted line shows rioo. 
right: The LF of Hydra A observed with 2MASS photoinetric data in the A'-band within rioo in absolute magnitude form, including background galaxies 
(black diamonds) and that the subtracted background contribution (red crosses). The vertical dashed line represents the complete limit of 2MASS. Blue 
and green solid lines show the best-fit Schechter functions when L* was fixed to the result from the stacked LF of galaxies in clusters by Lin, Mohr, 
Stanford (2003) and allowed to vary, respectively. 



to the virial radius in subsection 5.5 to constrain the star for- 
mation history in the cluster 

5.1. Comparison of temperature and entropy profiles with 
other systems 

Figure 12 compares the temperature profiles of the Hydra 
A cluster and several clusters observed with Suzaku, scaled 
with average temperature, < kT >, and ?'2oo- Here, r2oo = 
2.77/iTO yJ{kT)ll{) kcV /E{z) Mpc, with E{z) = (17^(1 + 
z)3 + 1 - r2ji/)i/2 (Henry et al. 2009). Then, raoo of the Hydra 
A cluster is 1.48 Mpc. Within 0.5r2oo, the temperature profiles 
agree with the scaled temperature relationship by Pratt et al. 



(2007), which is derived from a sample of 15 clusters observed 
with XMM out to ^ r^oQ. In contrast, beyond 0.57'200i the ob- 
served temperatures scatter by a factor of two, and tend to be 
lower than the XMM relationship. In the scaled temperature 
profiles, there was no evidence of any significant dependence 
on the average ICM temperature, or system mass. Some vari- 
ation in the temperature profiles should reflect the azimuthal 
variation as observed in Abell 1689 (Kawaharada et al. 2010) 
and Abell 1246 (Sato et al. 2012). 

In figure 13, we compare the entropy profile of the Hydra 
A cluster with other clusters observed with Suzaku; these pro- 
files are scaled with the average ICM temperature and r2oo by 
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Fig. 12. Radial profiles of temperature of clusters observed with 
Suzaku and XMM. Here, the temperature profiles are scaled with 
the average ICM temperature, < kT >, and are azimuthally aver- 
aged. We plotted the radial profiles of cluster of galaxies Abell 1689 
(Kawaharada et al. 2010), Abell 2142 (Akamatsu et al. 2011), Abell 
1413 (Hoshino et al. 2010), PKS-0745 (George et al. 2009), the Perseus 
cluster (Simionescu et al. 2011), Abell 1246 (Sato et al. 2012), the 
Virgo cluster (Urban et al. 2011), RXJ1159 (Humphrey et al. 2011) 
and the Hydra A cluster (current study). The scaling radius r'2oo was 
calculated by Henry et al. (2009). The solid line represents the scaled 
temperature profile by Pratt et al. (2007). 

Henry et al. (2009). Here, we used the weighted average of the 
radial profiles of different directions of each cluster Contrary 
to the expected r^-^ relationship, these profiles become flat be- 
yond 0.5r2oo. Except for PKS-0745 (George et al. 2009), the 
scaled entropy profiles of these clusters were universal with 
ICM temperatures above 3 keV and did not depend on the ICM 
temperature. For the PKS-0745 cluster, Eckert et al. (2011) 
determined that the ROSAT surface-brightness profile is sta- 
tistically inconsistent (7.7cr) with Suzaku results beyond r2oo 
and that the difference is likely explained by the existence of 
additional foreground components at the low Galactic latitude 
of the source. These components were not considered in the 
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Fig. 13. Radial profiles of entropy of cluster observed with Suzaku 
and scaled: with the average ICM temperature, < kT >, and ai'c az- 
imuthally averaged. Data for Abell 1689 (Kawaharada et al. 2010) are 
plotted in dark blue, for Abell 2142 (Akamatsu et al. 2011) in green, 
for Abell 1413 (Hoshino et al. 2010) in orange, for PKS-0745 (George 
et al. 2009) in puiple, for Perseus (Simionescu et al. 201 1) in light blue, 
for Abell 1246 (Satoet al. 2012) in red, and for Hydra A (current study) 
in black. The scaling radius T200 was calculated by Henry et al. (2009). 
The solid line shows the power-law profile with slope of 1 . 1 . The nor- 
malization of the power-law is scaled with our result for the Hydra A 
cluster. 

Suzaku background modeling. 

5.2. Ion-electron relaxation in the Hydra A cluster 

One interpretation for the flattening of the entropy profile 
is that caused by deviations in electron and ion temperatures 
(Hoshino et al. 2010; Akamatsu et al. 2011). The time-scale 
for thermal equilibration between electrons and ions through 
Coulomb scattering is given by 

l.SkcvJ 



i. 



'^■^^ilO-4cm-3J 



Gyr 



(8) 



(Takizawa & Mineshige 1998; Takizawa 1998; Takizawa 1999; 
Akahori & Yoshikawa 2010). Considering this time-scale and 
the shock propagation speed, the radial length above which the 
electron temperature is significantly lower than the ion tem- 
perature is proportional to the square of the ICM temperature 
(Takizawa 1999). On the other hand, r2oo is proportional to the 
square root of the ICM temperature. Therefore, some depen- 
dence on the mean ICM temperature is expected in the tem- 
perature and entropy profiles. However, there is no system- 
atic dependence on the average ICM temperature in the tem- 
perature and entropy profiles, as shown in figure 12 and figure 
13. Furthermore, if the ICM temperature is underestimated and 
the actual entropy profile of the Hydra A cluster follows the 
r^-^ relationship, a flat ICM temperature profile is necessary, 
i.e., r~" out to rioo, considering the electron density profile 
of J--!! 6~i 7) beyond 400 kpc. However, the flat temperature 
profile contradicts to the results of numerical hydrodynamical 
simulations, where negative gradients of the temperature pro- 
files outside cool core regions are naturally produced (Borgani 
et al. 2002; Burns et al. 2010; Nagai & Lau 2011). Therefore, 
for the Hydra A cluster, it is difficult to explain the flattening 
of the observed entropy profiles by deviations of ion-electron 
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temperatures. 

5.3. Deviations from hydrostatic equilibrium 

The remaining interpretation for the flattening of the en- 
tropy profiles at the cluster outskirt is deviations from hydro- 
static equilibrium, because when the gas-clumping effect is 
significant, we also expect that the ICM also deviates from 
hydrostatic equilibrium. Then, infalling matter may have re- 
tained some of its kinetic energy in the form of bulk motions, 
and the thermal energy deficit in the ICM yields lower en- 
tropy. Furthermore, the underestimate of the gravitational mass 
and/or gas clumping leads to an overestimation of the baryon 
fraction. Recent numerical simulations by Vazza et al. (2009) 
and Nagai, Vikhlinin, Kravtsov (2007) have shown that the ki- 
netic energy of bulk motion carries ^ 30% of the total energy 
around the virial radius. With simulations, Nagai & Lau (201 1) 
showed that beyond r2oo, gas clumping leads to an overestima- 
tion of the observed gas density and causes flattening of the 
entropy profile. If the hydrostatic mass was underestimated by 
130% in the filament direction and 100% in the void direction, 
the baryon fractions at 7'ioo would have the same value as that 
of WMAP7, although the observed gas fraction calculated from 
the hydrostatic mass in the void direction is consistent with the 
WMAP7 result within error bars. 

The distributions of the ICM and galaxies in the Hydra A 
cluster elongate in the filament direction. If the ICM is un- 
der hydrostatic equilibrium in a non-spherical dark matter halo, 
the X-ray emission and dark matter halo elongate in the same 
direction. Numerical simulations by Burns et al. (2010) de- 
termined that gas, galaxies, and dark matter continue in the 
filament direction due to the accretion of gas and subclusters. 
Then, it is reasonable that the dark matter halo also elongates 
in the filament direction. However, the hydrostatic mass in the 
void and filament directions of the Hydra A cluster are in good 
agreement. These results indicate that deviations from hydro- 
static equilibrium should be more significant in the filament 
direction than in the void direction. 

The clumping effect will also be higher in the filament di- 
rection, since smaller systems are thought to be accreted from 
this direction. For the Hydra A cluster, we were not able to 
detect clump candidates with Suzaku. The spectra of detected 
point sources are consistent with those of background active 
galactic nuclei. Considering the luminosities of detected point 
sources, the luminosity of clumps should be smaller than a few 
lO''^ erg s^^, which corresponds to cores of groups of galax- 
ies. Therefore, if clumps are significant, their scales should be 
smaller than those of small groups of galaxies. 

5.4. The concentration-mass relationship 

From numerical simulations, a weak variation in the concen- 
tration is expected from low-mass to high-mass clusters, re- 
flecting differences in the formation epochs of low-mass and 
high-mass halos (Navarro et al. 1997; Bullock et al. 2001; 
Dolag et al. 2004). Pointecouteau, Amaud, Pratt (2005), Ettori 
et al. (2011) and Okabe et al. (2010) investigated the relation- 
ship between the concentration parameter and the cluster mass, 
or the concentration-mass (c — M) relationship. The results 
are compatible with the intrinsic dispersion of theoretical pre- 
dictions. Figure 14 compares the c — M relationship of the 
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Fig. 14. Observed distribution of concentration parameters as a func- 
tion of cluster masses. The mass of weak lensing (blue) and X-ray (light 
blue) were defined by Okabe et al. (2010) and Ettori et al. (2010);Ettori 
et al. (2011), respectively. Okabe et al. (2010) used A/no and cno, 
while Ettori et al. (2010);Ettori et al. (201 1) used M200 and C2oo- We 
considered J\/ioo for Hydra A. The solid line indicates the median re- 
lationship determined through CDM simulations for the WMAP5 cos- 
mological model, while the region enclosed within the dashed fines cor- 
responds to tfie range of (T(logioc) = 0.1, witfiin wfiicfi most of sim- 
ulated clusters are distributed (Duffy et al. 2008). The large diamond 
shows the observed mass profile of Hydra A (current work) assuming 
hydrostatic equilibiium. Black and red open tiiangles coiTespond to 
cases in which the /gas(< t) value of Hydra A (current work) is the 
same as the cosmic mean baryon fraction (Komatsu et al. 2011) in the 
filament and the void directions, respectively. 

Hydra A cluster and other clusters observed with weak-lensing 
(Okabe et al. 2010), and clusters observed with XMM by Ettori 
et al. (201 1). An assumption of cluster dynamical state or hy- 
drostatic equilibrium is not required in the weak-lensing re- 
sults. The distribution of the c — M relationship of these clus- 
ters agrees well with the relation expected from the numerical 
simulations by Duffy et al. (2008), although Ettori et al. (201 1) 
found that relatively low-mass systems tend to have higher con- 
centration parameters. 

The hydrostatic mass in the Hydra A cluster of the filament 
and void directions are represented with the same NEW mass 
model with Cvir = 9-10 and Afvir = (2.2 ± 0.4) x lO" Mq. 
As shown in Figure 14, this Cvir — Afvir relationship of the 
Hydra A cluster was outside a{\ogn^c) = 0.1 within which 
most clusters from numerical simulations are distributed Duffy 
et al. (2008). If the hydrostatic mass was underestimated and 
the baryon fraction using the gravitational mass beyond r,5oo 
of the Hydra A cluster is the same as the cosmic mean baryon 
fraction in the two observed directions, the fits with the NEW 
model yielded smaller values of Cvii and higher values My-„ as 
shown in Figure 14. In this case, the derived value of A/vii in 
the filament direction became higher than that in the void direc- 
tion. This result also indicates that the ICM in the outskirts of 
Hydra A deviates from hydrostatic equilibrium, and the X-ray 
emission and dark matter halo elongate in the same direction. 

5.5. Metal-Mass-to-Light Ratios out to the virial radius 

We derived integrated IMLR, mass-to-light ratio for Fe, out 
to the virial radius for the first time. Here, we show the IMLR 



14 



T. Sato et al. 



[Vol. 



Ooma 

Centaurus 

AWM7 

Hydra-A 

A262 



NGC5044 
NGC1550 




Fig. 15. Compaiison of iron-mass-to-light-ratio for other clusters. The 
radius is scaled with the characteristic radius r2()0 (Henry et al. 2009). 
For comparison, we plotted the radial profile of cluster of galaxies 
AWM 7 (Sato et al. 2008, Sato et al. in preparation), Abell 262 
(Sato et al. 2009b), the Centaurus cluster (Sakuma et al. 201 1) and the 
Coma cluster (Matsushita et al. 2011) and galaxy groups NGC 5044 
(Komiyama et al. 2009) and NGC 1550 (Sato et al. 2010). Open cir- 
cles are the result from XMM-Newton satellite, and others are Suzaku 
satellite. 

profile for the whole azimuthal angle to compare other clusters 
observed with Suzaku. Within 0.3 r2oo and beyond 0.3 r2oo, 
we used the Fe abundances obtained by Matsushita (20 11) and 
this work, respectively. Here, within 8'-15', or 0.35-0.65 r2oo, 
we used the weighted average of the Fe abundance at the three 
radial bins in this region. The error bars of the mass-to-light 
ratio include the errors in both the abundance and the K-band 
luminosity caused by the uncertainty in the contribution of 
fainter galaxies below the 2MASS threshold (see section 5.5 
for details). Figure 15 shows the cumulative IMLR profile of 
the Hydra A cluster The IMLR profile increases with radius 
within 0.4-0.5 r2oo, and becomes flatter beyond 0.5 7*200 with 
the value of ^ 0.01 Mq/Lk,q, which reflects the flattening of 
the gas-mass-to-light ratio beyond r5oo. For the void and fila- 
ment directions, the IMLR profiles should increase with radius 
beyond r5oo. 

At r2oo, the cumulative IMLR was 0.007-0.01 Mq/Lk,q- 
A previously reported theoretical model predicts that the 
oxygen-mass-to-light ratio (OMLR) of a cluster is a sensi- 
tive function of the slope of the initial mass function (IMF) 
(Renzini 2005). Here, the oxygen mass is a sum of that 
trapped in stars and that in the ICM. By adopting a Salpeter 
IMF with a slope of 2.35, the expected value of the OMLR 
is ~ 0.1 Mq/Lb,q- In contrast, a top-heavy IMF with a 
slope of 1.35 overproduces metals more than that with a fac- 
tor of 20. For the Hydra A cluster, by adopting the ra- 
tio of stellar-mass-to-light ratio 4-6 in the B-band and K- 
band (Nagino & Matsushita 2009; Matsushita et al. 201 1), the 
IMLR in the B-band becomes 0.03-0.06 Mq/Lb,q and the 
cumulative OMLR within the virial radius becomes ^0.1-0.3 
[Aq / A-pc) Mq / L B ,Q- Here, Aq and A-p^ are abundances of O 
and Fe in solar units, respectively. Beyond rsoo of the Hydra 
A cluster, the integrated stellar mass is only several percent of 
the integrated gas mass. Thus, even by considering the differ- 



ence of stellar metallicity and the abundance of the ICM, most 
of the metals in the cluster are in the ICM. Therefore, for cases 
with solar O/Fe ratio, the OMLR is similar to the expectations 
reported by Renzini (2005), assuming a Salpeter IMF in the 
Hydra A cluster In contrast, by adopting the O/Fe ratio from 
a nucleosynthesis model of SN II with metallicity = 0.02, as 
reported by Nomoto et al. (2006), the OMLR within the virial 
radius becomesO. 3-1. Oi\/0/iB.0 in the Hydra A cluster This 
value is larger than the expectation from the Salpeter IMF, and 
is consistent with a flatter IMF slope. 

In Figure 15, the cumulative IMLR profile of the Hydra A 
cluster was compared with those of other groups and clusters, 
including Coma (8 keV; Matsushita et al. 2011), Centaurus 
(4 keV; Sakuma et al. 2011), and AWM 7 clusters (3.6 keV; 
Sato et al. 2008), as well as NGC 1550 (1.2 keV; Sato et al. 
2010) and NGC 5044 groups (1.0 keV; Komiyama et al. 2009). 
Within 0.5r2oo, the cumulative IMLR profile of the Hydra A 
cluster increased with radius by a factor of two higher than for 
other systems. No systematic dependence on the ICM tem- 
perature was evident for clusters with temperatures higher than 
2 keV, although beyond 0.1r2oo, a group of galaxies shows a 
significantly smaller IMLR 

The stellar and gas-mass fractions within r5oo depend on the 
total system mass (Lin, Mohr, Stanford 2003; Lin & Mohr 
2004; Vikhlinin et al. 2006; Sun et al. 2009; Giodini et al. 
2009). These studies determined that within r5oo the stellar- 
to-total-mass ratios of the groups are much larger than those in 
the clusters, whereas the gas-mass fraction increases with the 
system mass. The observed higher stellar mass fraction and the 
lower gas mass fraction within rsoo in poor systems are occa- 
sionally interpreted as proof that the star formation efficiency 
depends on the system mass. However, as shown in Figure 
11, the gas is more extended than the stars in a cluster within 
rsoo. The gas density profiles in the central regions of groups 
and poor clusters were observed to be shallower than those in 
the self-similar model, and the relative entropy level was corre- 
spondingly higher than that in rich clusters (Ponman, Cannon, 
Navarro 1999; Ponman, Sanderson, Finoguenov 2003; Sun 
et al. 2009). Then, the difference in the ratio of gas-mass-to- 
stellar-mass may reflect differences in distributions of gas and 
stars, which in turn reflects the history of energy injection from 
galaxies to the ICM. To study the fractions of stars and gas in 
clusters of galaxies, and its dependence on the system mass, 
we need measurements of gas and stellar mass beyond rsoo of 
other clusters. 

6. Summary and Conclusion 

We presented the results of the Suzaku observation of the 
Hydra A galaxy cluster, which is the first medium-sized clus- 
ter (temperature ^ 3 keV) observed with Suzaku out to the 
virial radius. Two observations were conducted, north-west 
and north-east offsets, which continue in a filament direction 
and a void direction of the large-scale structure of the Universe, 
respectively. We investigated possible azimuthal variations 
in temperature, electron density, entropy and total mass pro- 
files. Our analysis revealed that distributions of X-ray emis- 
sion and galaxies elongate in the filament direction. The en- 
tropy profiles become flatter beyond rsoo in contrast to the 
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r^^^ relationship expected from the shock heating model of 
the ICM. The entropy profiles are universal in clusters ob- 
served with Suzaku when scaled with the averaged tempera- 
ture of each cluster The hydrostatic masses in the two direc- 
tions agree and the NFW universal matter profile represents 
the hydrostatic mass distribution obtained up to 1.8 r5oo, with 
Cvir = 9-10 and Afvir = (2.2 ± 0.4) x IO^^Mq. The gas frac- 
tion, Afgas(< ?')/AfH.E.(< r) significantly exceeds the cosmic 
mean baryon fraction of WMAP7 beyond ^ rsoo- The flat- 
tening of the entropy profile and higher gas fraction contra- 
dict expectations based on numerical simulations. An under- 
estimate of gas temperature due to the discrepancy between 
ion and electron temperatures gives higher entropy and smaller 
gas fraction at cluster outskirts. However, a flat or increasing 
temperature profile is required. If bulk motions caused by in- 
falling matter from filaments of clusters are higher than that of 
numerical simulations, and if the ICM deviates from the hy- 
drostatic equilibrium, the entropy becomes smaller and the gas 
fraction is overestimated. In addition, we obtained IMLR up 
to the virial radius for the first time, and compared our results 
with other clusters and groups. The IMLR profile is consistent 
within 0.5 r2oo with other clusters, and become flatter from 0.5 
r2oo tor2oo- 
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